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TECHNICAL MEMORANDUM 


PRELIMINARY ANALYSIS TECHNIQUES FOR RING AND 
STRINGER STIFFENED CYLINDRICAL SHELLS 

I. INTRODUCTION 


Over the years many methods of stiffened panel analysis and design have been proposed and 
substantiated by test. Most methods deal with only one element of the stiffened structure, such as 
the skin, the stringer, or the ring. Not many references combine all the elements of skin stringer 
design and analysis into one cohesive process. This is especially true when the skin of a structure is 
allowed to buckle prior to application of the ultimate load. This report outlines methods of analysis 
for the major failure modes for the buckling of thin-walled circumferentially and longitudinally 
stiffened cylindrical shells. The report is intended particularly to address launch vehicle design 
issues. Loading on the vehicle will consist of pure bending, axial compression, and shear, all in the 
elastic range. Generally, any advance of the load beyond the buckling limit is considered a structural 
failure and must be avoided in launch vehicle design. (The skin, however, may be allowed to buckle 
at limit loads.) A Microsoft Excel worksheet with accompanying macros has been developed to 
facilitate application of the various analysis methods. These analysis programs are available by 
request from the author. 

The analysis methods presented are organized according to failure mode. All necessary 
design curves have been curve fit to allow automated analysis in the spreadsheet program. Sections 
detailing the calculation of stress in the vehicle, as well as calculation of margins of safety, are also 
included in the paper. The appendices contain hand calculations, additional analysis information, and 
the analysis programs. This report will focus on the integrally Tee stiffened shell. The reader should 
note that, unless otherwise stated, all methods presented in this report are for use in the elastic 
region. 


II. ANALYSIS METHODS 


Analysis of the stiffened shell begins with recognition of the various failure modes. The failure 
modes listed below encompass the most significant failure modes of the shell. Analysis techniques 
for each of the failure modes listed will be presented. 

Buckling failure modes can take one or more of the following forms. 

( 1 ) Classical bifurcation buckling 

(2) General instability 


(3) Stringers 

(a) Local buckling 

(b) Crippling 

(c) Column failure 

(4) Skin 

(a) Compression buckling 

(b) Shear buckling 

(c) Pressure. 

Frames, of course, may experience cap failures, web failures, and other buckling failures. 
However, the interest in frames at this level of design is to determine an acceptable moment of 
inertia, or other overall general characteristic of the frame design, which will stabilize the structure 
against general instability failure. For these reasons, detailed frame design is excluded from this 
report. Also, the primary focus of this report is application to pressurized shells. Pressurization pre- 
cludes penetration of the tank by rivets except in extreme cases. Therefore, inter-rivet failure and 
face sheet wrinkling have been excluded as failure modes. 


A. Bifurcation Buckling 

The difference between bifurcation buckling (also commonly referred to as “classical” or 
“classical bifurcation” buckling) and other types (or modes) of buckling failure are often confusing. 
The point (load value) at which a column fails due to bifurcation buckling represents the intersection 
of two equilibrium paths in the structure. The failure modes presented in subsequent sections of this 
report represent collapse, or failure, at a limit point. Figures la and lb illustrate the difference 
between bifurcation buckling and failure at a limit point. The variable P represents the applied load, 
and A represents the displacement. 




Figure 1. Bifurcation buckling. 
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Figure 1 a represents the load-displacement curve for a Bellville spring. The point on the 
equilibrium path at which load P is a relative maximum is called a limit point. A limit point could be 
column buckling, skin buckling, stringer crippling, etc. Figure 1 b represents the load-displacement 
curve for a thin-walled cylindrical panel under axial compression. In this figure, the “primary (or fun- 
damental) equilibrium path is intersected by a secondary path.” 1 The point of intersection is called 
the bifurcation point. 

Failure of a general shell usually occurs through collapse at some limit point rather than 
through bifurcation. However, bifurcation buckling of the shell must be considered in the design pro- 
cess. It is important that the analyst realize that “the classical (or bifurcation) buckling analysis 
may give results of little or no value if the shell geometry deteriorates appreciably (Brazier effect) or 
stresses are redistributed ... in the subcritical load range.” 3 In the examples used in this report, the 
skin of the cylindrical shell is allowed to buckle at the limit load while the stringers and effective skin 
are allowed to buckle only at ultimate load. There is considerable redistribution of stress. Therefore, 
the classical bifurcation solution alone is of little use, but will be calculated as an illustration of the 
method. 

Computation of the linear bifurcation buckling load and application of an empirical knock-down 
factor provides a conservative method of determining an appropriate allowable load level. The 
following analysis techniques combine the wide column allowables with the bifurcation buckling 
allowable reduced by a “knock-down factor” which is a function of (R/t) e . 

The bifurcation analysis does yield good results for bending of cylinders with (R/t) e values 
large enough so that the Brazier effect is negligible. It has been shown “that the use of wide column 
load as a design limit for stringer-stiffened cylinders was unduly conservative. It was suggested that 
a term be added to the wide-column load which corresponds to the curvature effect. This term was 
obtained as the difference between the classical buckling load and the wide column load multiplied by 
a reduction factor.” 3 

The effect of curvature is introduced by taking the difference between the wide column and 
classical allowables, multiplying by a “knock-down factor” ( <p ), and adding this result to the wide 
column allowable. 3 


Ncr — A'Wc + <p (N CL -N WC ) 


( 1 ) 


1. Classical 

Presented below are the constitutive relationships for the orthotropic shell. 1 Examples of 
construction that may be treated as orthotropic include corrugated sheets, fiber reinforced plastic 
sheets, and plates with closely spaced stiffeners. This method yields adequate results for closely 
spaced rings. For spacing greater than 30 to 40 in, setting the ring properties equal to zero may yield 
an adequate solution. 1 
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where [C,y] represents the stiffness matrix for the shell. Definition of the stiffness parameters is 
given below: 4 
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where 


J = torsional stiffness constant 

A = area of stringer or ring 

/ = area moment of inertia of stringer or ring 

e = distance from the skin middle surface to the centroid of the stiffener cross section 
b S k = stringer spacing. 

The coupling parameters C14 and C25 are positive for stiffeners outside the skin and negative 
for stiffeners inside the skin. 
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The buckling coefficients are written in terms of the axial half wavelength number (m), and 
the full circumferential wave number (n): 

A(m)=[m ff = A„ , = 


From reference 4, the matrices oq and 31-6 defined as. 


a 0 (m,n ) = 


A n Aj 2 

a 2] a 22 


and 


where. 


a x (m/i) 


A u A 12 A 13 

■^12 ^22 ^23 > 

A 13 A 23 A 33 


A n (m,n) = C n A m +C 33 £ „ , 
A n (m,n) = [C 12 +C 33 ]m ( £■) ^ * 
A 22 (m,n) = C 33 A m +C 22 £ „ , 


4 13 (m,n) = m (f ) + c ij^ (-^) + [C is+JCjJm (^) £ „ , 

A 23 (m,„) = [C I5+ 2C 36 ]A „(#)+% (f ) . ($) , 

■M»vO = C„A J,4{C« + 2C 45 ]A „+C 55 | J+^+2(^ „+2(%)A . . 


Also from reference 4, the basic buckling equation is defined as, 


m + Ny% n 


|a 1 (m,n)| 
\a 0 (m,n)\ ' 


Here, the effect of internal pressure is included by calculation of the pressure induced line load in 
pounds per inch ( N y = - p*R ), where p is internal gauge pressure. N x is then determined by the equa- 
tion 


N x = 



aM,n)\ ' 

\a 0 (m,n)\ y " 


The classical bifurcation buckling load ( Ncl ) is determined by attempting all combinations of m and 
n. Ncl is the minimum of the buckling values obtained. 
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2 . 


The wide column buckling allowable is obtained from the matrix equations by the following 
relation : 3 

A/ tvc = /? 2 C 4 4^j , 

C 44 = C 44 — and C 55 = £ 55 — 7 ^ • 

'-'11 '-'22 

In order to combine the wide column and classical buckling solutions, the “knock-down fac- 
tor” ( 9 ) is required: 

<P=fW) e . 

The “knock-down” factor (<j£>) as a function of (R/t) e is obtained from figure 2 . 3 A probability 
level of 99 percent is recommended. 


1.0 


as 

0.2 


0.1 

10 20 SO too 200 SOO 10C0 2000 

(R/t > e 

Figure 2. Empirical “knock-down” factors. 

W) e = 


All quantities are now known and can be applied to the relation: 

Ncr = N w c*-(p(N C i-N wc ) . 
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B. General Instability 


The purpose of general instability calculations is to avoid general instability failure, as illus- 
trated in figure 3. 1 * * * 5 Frames are designed to preclude general instability failure, or rather to ensure 
panel failure as illustrated in figure 4. 5 


A 



of Bending Spring Behavior of Frames 


Figure 3. General instability buckling. 



of Bending Spring Behavior of Frames 

Figure 4. Panel instability buckling. 


1. Shanlev Criteria for Cylindrical Shells in Bending 

To prevent general instability, Shanley has determined an expression— equation (3)— for the 
required product of frame modulus of elasticity and moment of inertia for pure bending of a stiffened 

shell: 5 


{ET) f = C f MD 2 /L . (3) 

The coefficient C/has been determined through experimentation to be 62.5xl0' 6 (or 1/16,000). 

Figure 5 shows the data from which this coefficient comes. One can observe the crossover point from 
general instability failures to panel failures when the value of C/is approximately equal to 62.5x1 O' 6 . 
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160 


Cf xio 6 

Figure 5. Frame buckling coefficient. 

Of course, not all loading in launch vehicles is pure bending. This problem may be remedied by 
calculation of an equivalent moment where, 6 

M eq = PR/2 , (4) 

where P represents the axial load on the cylindrical shell. The equivalent moment is combined with 
the pure moment to get total effective moment which is then returned to equation (3) for calculation 
of the required frame (El). 

It should be pointed out that the Shanley method of frame sizing may give either conservative 
or unconservative results. The results depend on the configuration under analysis. The Shanley 
method should only be used as an initial sizing measure. Correlation of results with other methods 
would be advisable. 
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where 


b = stringer spacing 

d = frame spacing 

R = cylinder radius 

t = skin thickness 

A SI = stringer area 

Af = frame section area 

t s = distributed stringer area = A st /b 

tf = distributed frame area = A/d 

If = bending moment of inertia of frame section 

If = distributed frame moment of inertia = If/d 

p s = stringer section radius of gyration 

Pf = frame section radius of gyration 

L = length of cylinder 

E = modulus of elasticity. 

If the frames are not attached to the skin, the coefficient 4.8 in equation (6) should be 
replaced by 3.25. The effective skin width for frames should be taken as the total frame spacing. 5 
The effective width for stringers for use in frame calculations is given by the following equation. 5 
Effective skin is that skin which is assumed to act with an adjoining element and carries the same 
stress as that element. 


w (Fc \' k 

^=0' 5 bf) ■ 5) (7) 

Fc cr is the critical buckling stress for a curved skin panel, and F c is the compressive stress at bend- 
ing general instability — applied ultimate stress. 

A distinct advantage of using the Shanley method is that design parameters such as stringer 
geometry and skin thickness do not have to be known. Use of the Becker method, however, requires 
that a preliminary design exist for evaluation including stringer and skin definition. Also, there are 
several ambiguities in the Becker equation that Dr. Bruhn does not clarify. These ambiguities stem 
from the definition of sectional properties, and whether or not to include effective areas and the like. 
After review of the original reference, an equivalent equation for critical stress can be obtained which 
more clearly defines the use of effective skin for the stringers and frame. 

° cr~ CE c Q b , (ref. 7) (8) 


where 


C = 4.800 for frames attached to the skin 
C = 3.25 for frames not attached to the skin 


9 


Ub- R 

In determining the radius of gyration for the stringers and frames, the effective widths are 
used. The effective width for the frames is defined as the frame spacing itself, while effective width to 
be used with the stringers is calculated by equation (7). According to Becker, these equations are 
valid only in the range for: , 


Good agreement can be had between the Becker and Shanley methods even for values of 
L 2 /RT below 100. The frame section properties, however, should not be increased by the standoff 
distance of the frame from the skin. The properties should be calculated as if the frame were adjacent 
to the effective skin. For critical values of stress in the plasticity region, the secant modulus can be 
substituted for Young's modulus with credible results. 


C. Stringers 


1. Local Elastic Buckling 

“Thin flat sheet is inefficient for carrying compressive loads because the buckling stresses 
are relatively low. However, this weakness, or fault, can be greatly improved by forming the flat 
sheet into composite shapes such as angles, channels, zees, etc.” 5 Calculation of the composite 
buckling strength is thus necessary to prevent failure of the stringer column. However, since the 
stringer will continue to carry load after local buckling has occurred, local buckling may be allowed in 
some instances. Local buckling is more likely to be a design driver where substantial deformation of 
the stringer flanges causes debonding of insulation or other material. In these instances, it is accept- 
able to compare the local buckling allowable to the limit load for margin of safety calculations. 

Analysis of the local elastic buckling failure mode is easily accomplished by dividing the 
flanges of the section into individual plate elements with large a/b ratios. Jumping ahead to buckling 
of the skin, and using the buckling equation for a flat plate, equation (26), 

_ k c * 2 E I,y 
c 12(l-u 2 )VfcJ ' 

Figure 1 1 is used to determine the buckling coefficient k c . (Figure 1 1 and the buckling equa- 
tion are explained in detail in the skin buckling section of this report.) Since at least two of the 
flanges are usually of equal size, they buckle at the same stress. Therefore, they cannot be relied 
upon for edge support. For this reason, a simply supported edge condition is assumed along the 
longitudinal junction of the flange elements. If the opposite edge is free, the buckling coefficient is 
0.43. If the opposite edge is also simply supported, choose k c equal to 4.0. 

The flange width b extends to the centerline of the adjacent leg for formed angles. For 
extruded angles, the width b extends to the inside edge of the adjacent flange or leg. 5 The smallest 
buckling stress found in the composite shape, not the average, becomes the critical buckling stress 
for local elastic buckling. 


10 



2. Crippling 


“Tests of short lengths of sections composed of flange-plate elements often show that after 
the section has buckled locally, the unit still has the ability to carry greater loads before failure 
occurs. . . . For cases where local buckling occurs at low stress, the crippling or failing stress will be 
higher. When local buckling occurs at high stress, such as 0.7 to 0.8 Fey, buckling and crippling 
stress are practically the same.” 5 For clarification, stringer crippling may be viewed as a material 
failure where critical values are compared against material ultimate or yield stress. Local elastic 
buckling is considered a stability failure. 

Three methods for computation of the stringer crippling allowable are presented here. The 
method from the NASA Structures Manual is the simplest of the methods to employ. 8 


a. NASA Structures Manual . The NASA Structures Manual provides a detailed step-by- 
step procedure for determining the overall strength of a sheet and stiffener combination. The method 
is very similar to the Gerard method, but with some modifications. No definite reference for this 
method was given. However, all references listed by the structures manual were dated earlier than 
the paper delivered by Gerard which details his method of stringer crippling determination. 

The stringer crippling stress is determined by the following equation: 8 


F = 

1 CS 


Zbjn 


(9) 


One can easily see that this stress represents the average failing stress of the stringer ele- 
ments or flanges. The failing stress of individual elements is f cc „. Dimensions of the individual ele- 
ments are determined consistent with figure 6. 8 Note that two number 1 elements will be needed for 
analysis of the Tee stringer. The failing stress of each individual element is found from figure 7. 8 


h^-i 



n 





Figure 7. Nondimensional crippling curves. 










Following the procedure outlined in the manual, the effective width is not used in determina- 
tion of the crippling stress of the stringer. This seems contrary to assumptions of the stress distri- 
bution in the shell because at this point the effective skin acts with the stringer; experiencing the 
same stress and load. Common practice is to include the effective width as another flange element of 
the stiffener when performing crippling stress calculations. The example problem will include the 
effective skin as part of the stringer in determining crippling stress. 

b. Needham . The Needham method is most useful for formed or extruded stringers that are 
mechanically attached to the skin — such as a hat stringer riveted to the skin. The method consists of 
dividing a stringer into angle sections. The strength of each of the angle sections is determined, and 
then the total strength of the stringer is achieved by summing the individual section strengths. 
Needham has arrived at equation (10) for determination of angled section strengths: 5 


where 


F C 

* Cs 


( 10 ) 


b'/tsk = equivalent b/t of section = ( a+b)/2t S k 

(a and b are the leg elements of the angle) 

C e = coefficient that depends on the edge support 
0.316 (two edges free) 

0.342 (one edge free) 

0.366 (no edge free). 

The crippling load of the angle may then be determined as follows: 

P CS = F CS A , (ref. 5) (11) 

where A is the area of the element in question. 

The total crippling stress of the stringer representing the average of all stringer angles is 

then: 


1 Crippling Loads of Angles 

*cs ~ v a e a i ‘ (ref. 5) (12) 

L Area of Angles 

c. Gerard . The Gerard method can be thought of as a broader application of Needham. The 
crippling stress equations for various stringer configurations are presented here. Equation (13) is for 
sections with distorted unloaded edges such as angles, tubes, V groove plates, multicorner sections, 
and stiffened plates. The accuracy of this equation is said to be ±10 percent, as reported by Gerard. 9 


= 0.56 [{gt 2 /Ap/F cy ) 1/2 ] 0 ' 85 . (ref. 5) 


(13) 


Equation (14) is required for sections with straight unloaded edges such as plates, Tee, 
cruciform, and H sections. Reported accuracy is within ±5 percent. 
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(ref. 5) 


(14) 


^ = 0.67[(«( J MXHF ( ,) ,n ] OJ . 

For two comer sections, J, Z, and channel sections, use equation (15). Accuracy is within ±10 
percent. 

^ = 3.2 [{t 2 /AXE/F cy Y n f - 75 . (ref. 5) (15) 

Equations (13) through (15) represent approximations or simplifications of the data pre- 
sented by Gerard. The general equation for stringer crippling is given by Gerard as: 

Tr = P,[(slJ«MlEIF e) ) m ] m . (ref. 9) (16) 

The coefficient fi g is determined experimentally and results are tabulated by Gerard as func- 
tions of 7 Jt s k- Using P g - 0.56 would be considered an average for stiffened plates. The actual values 
range from 0.562 to 0.464, as shown in table 1. The constant g is the sum of the number of flanges of 
the angled elements and the number of cuts required to divide the stringer into angled elements. The 
exponent m is 0.85. 


Table 1. Gerard coefficients. 


^ w^sk 

K 

1.16 

0.562 

0.732 

0.505 

0.464 

0.478 


1 w is the average thickness of the flange sections as determined by equation (17): 


_ 


(ref. 8) 


(17) 


The skin thickness is represented by t S k, and the exponent m is determined experimentally. 
However, it does not change with ljt s k, and is dependent on the type of stringer arrangement. 


The data presented above were for T-stiffened panels. Inspection of the integrally stiffened 
panel indicates that it closely duplicates the Z-stiffened panel. Coefficients for Z-stiffened panels 
are: m = 0.85, and fi g = 0.558. fi g data as a function of ljt s k are the same as T-stiffened panels. The 
number of flanges and cuts will change. 

Bruhn also presents a series of illustrations in which both methods are used to determine 
stringer crippling stress. Depending on the stringer configuration, there can be significant differences 
in the Needham and Gerard methods. Also, the crippling value is subject to upper limits that should 
not be exceeded unless test data can substantiate such a move. A table of upper limits is presented 
as table 2. 5 
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Table 2. Maximum crippling stress. 


Type of Section 

Maximum F cs 

Angles 

0.7 F cy 

V Groove Plates 

Fcy 

Multicomer Section, Including Tubes 

0.8 Fcy 

Stiffened Panels 

Fcy 

Tee, Cruciform, and H Sections 

0.8 F cy 

Two Comer Sections, Zee, J Channels 

0.9 Fcy 


The reader should note that the entire width of skin between stringers is used as a flange in 
determining the stringer crippling value (using Gerard's method), as opposed to using the effective 
width only. This convention is maintained since the empirical equations were obtained using that 
convention. 

3. Column Failure 

In general, column failure is the limiting failure mode for most longitudinally stiffened vehicle 
structures. The primary buckling equation for elastic failure is simply stated as equation (18); the 
Euler buckling equation: 5 


7 _ k 2 E 
‘ (Up ) 1 ' 


(18) 


Equation (18) can be rewritten involving the tangent modulus for stresses in the inelastic region. 
Where E, is approximately E for stresses in the elastic range: 


F = 


k 2 E , 
(Up ) 2 


(19) 


The radius of gyration (p) for the column (or stringer), is calculated by equation (20): 5 

p = VT/A . (20) 


The buckling strength of a column is also heavily influenced by the end restraint on the 
column. Adding the end-fixity coefficient c into equation (19) allows incorporation of the end-fixity 
constraint into the buckling equation. A new effective column length is determined by equation (21): 5 

L' = {UJc) . ( 21 ) 


The Euler equation thus becomes: 


C C L'/p ) 2 ' 


( 22 ) 
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Determination of the tangent modulus is accomplished through use of the basic Ramberg- 
Osgood relationship. 5 Note that E = E, for stresses in the low elastic range. 



n = 1 + 


In ( 17 A) 

ln MaJ ’ 


( 23 ) 


where O0.7 is the secant yield stress found by drawing a line on the material stress-strain curve from 
the origin with a slope of 0.7E and reading the stress at intersection with the stress-strain curve. 
Ob.85 is found similarly. 

End-fixity coefficients for various end constraints are presented as figure 8. These coeffi- 
cients are presented for completeness. In practice, c = 1 or 1.5 is generally applied to the skin- 
stringer design problem. A value of 1.0 would be conservative. 



C=1 C=4 C=9 C=4 C=l/4 C=2. 05 

Figure 8. End-fixity coefficients. 


The NASA Structures Manual presents design charts (appendix D) for determining the end- 
fixity coefficients based on the bending stiffness of the end restraint. This data can be used in 
determining the slenderness ratio (L'/p) of the column. However, the analyst may be inclined to 
determine his or her own coefficient depending on the fidelity and conservatism desired of the analy- 
sis. 


Figure 9 shows a typical curve of F cr (critical buckling stress) as a function of L'/p. Buckling 
of columns with stable cross sections such as tubes and other closed sections follows the curve 
ABFC. Equation (22) is applied to determine the critical buckling stress. 

However, for columns of unstable cross section such as channels, Tees, and Tee stringers 
with effective skin width, curve DEFC must be employed. Critical stresses in the FC region are 
found through the Euler equation — equation (22). If the slenderness ratio of the column shows it to 
be in the transition region of the curve of figure 9, the buckling limit will be below that predicted by 
simple Euler buckling. 
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Figure 9. Typical column buckling curve. 


Stresses in the transition region are determined using the Johnson-Euler equation (24), 
where F cs is the crippling strength of the stringer: 5 


F = F =- cs — 
c a An 2 E 


(■r 


(24) 


To find out in which section of the buckling curve a column lies, it is necessary to compute the 
intersection point of the Euler and Johnson-Euler curves. Setting the buckling stresses from 
Johnson-Euler and Euler equations equal to one another and solving for L'/p, the intersection point is 
shown to lie at: 


(L'/p)=n<J 1 ^ . (ref. 5) (25) 

For values of L'/p greater than the intersection value, use the Euler equation. For values of L'/p less 
than or equal to the intersection value, use the Johnson-Euler equation. 

Now that methods for calculating the buckling strength of a column have been presented, it is 
necessary to apply these methods to the skin-stringer design problem. “As load is increased, the 
sheet buckles between the stiffeners and does not carry greater stress than the buckling stress for 
the skin. However, as the stiffeners are approached, the skin being stabilized by the stiffener to 
which it is attached can take a higher stress, and immediately over the stiffener the skin carries the 
same stress as the ultimate strength of the stiffener, assuming the sheet has a continuous connec- 
tion to the stiffener.” 5 
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In general, attempts to determine the amount of skin acting with the stringer have resulted in 
long and complex equations. To simplify the determination of the effective skin width, equations have 
been developed to determine the skin width that would be considered as taking a uniform stress 
equivalent to the stress in the effective stringer under actual nonuniform conditions. Figure 10 illus- 
trates the progression of panel buckling and the assumption of uniform stress made over the effective 
skin. 5 


A method for effective sheet width determination is presented here. The procedure follows 
the analysis of a flat sheet in compression with long edges simply supported. The critical buckling 
stress for the flat sheet is found from equation (26): 5 


F= k c* 2E (±) 2 

c I2(l-v 2 ) ybl ' 


(26) 


The k c value determined from experiment approaches 4.0 if the long edges of the sheet are assumed 
to be simply supported. 


Fig. b Sheet Stress Distribution Before Buckling 
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Figure 10. Progression of panel buckling. 

The problem of determining effective width has been simplified somewhat by Von-Karman 
and Sechlcr. 5 Their method consists of solving for an effective width ( w e ) in place of b in equation 
(26) with the critical stress equal to the yield stress of the material. However, since the buckling 
stress can be greater than the yield stress, the yield stress has subsequently been replaced by the 
stress in the stringer itself ( F s( ). Substituting 4.0 for k c and 0.3 for Poisson’s ratio, equation (26) 
reduces to: 
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or 


F c = 3.60 E(t/w e ) 2 
w e = l.90t(E/F st ) 05 - 


(ref. 5) 


(27) 


The “r” of equation (27) is defined as the thickness of the skin plus the skin and stringer land 
divided by two (see appendix A). 5 However, experiments by Newell indicate the constant 1.9 is too 
high, and 1.7 would be more appropriate, 5 possibly due to the conservatism in the buckling coeffi- 
cient. Although the preceding derivation was accomplished using a buckling coefficient of 4.0, a 
simply supported flat plate, the determination of effective width for curved plates does not differ 
significantly from flat plates for values of Z< 30. 10 

Z 4 (!-„’)*. 

The stringer should now be treated as a wide column made up of the stringer and its effective skin, 
and the appropriate Johnson-Euler or Euler equation applied to determine buckling allowables. 


D. Skin 


1. Compression Buckling 

The first evaluation of the skin strength comes in the form of simple flat sheet buckling. By 
using equation (26) and setting b to the stringer spacing and t equal to the skin thickness, the buck- 
ling strength of the skin between stringers can be determined: 

k^E , j _\2 

c 12(1 -d 2 )U/ ' 

Just as in determination of the effective width, the buckling strength of the skin is very much 
dependent on the buckling coefficient. A value of 4.0 is commonly used and is considered conserva- 
tive because of its assumption that the long edges of the sheet are simply supported. Figure 1 1 
illustrates the effect of increasing the a/b ratio on the k c values. 12 For long simply supported sheets 
(constraint C), the buckling coefficient value approaches 4.0. Figure 12 offers an alternative to this 
conservatism by allowing the engineer to choose a buckling coefficient that lies somewhere between 
the conservative simply supported case and the nonconservative clamped edge case. 13 

One should immediately notice that the buckling coefficient curves discussed thus far are for 
flat sheets, and the problem being investigated is one of stiffened cylindrical structures or curved 
sheets. Knowing this, however, many analysts use the flat sheet data instead of that for curved 
sheets, which are generally stiffer than their flat counterparts. For curved sheet panels, the buckling 
equation remains the same, however, k c is determined from figure 13. Curved sheets of large radius 
( b 2 /Rt < 1) can be analyzed as flat plates. 14 The data used in obtaining figure 13 are for a simply 
supported edge condition. 5 When curved sheet k c values are compared to those from figure 12, the 
flat sheet may be larger. This is caused by the simply supported edge restraint used to obtain figure 
13 data. It is recommended that the maximum k c resulting from figures 11, 12, and 13 be used in 
determining skin buckling coefficients. Data from figure 1 1 or k c = 4.0 may be used for conservatism. 
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Figure 12. Compression buckling coefficient. 
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Figure 13. Compression buckling coefficient for curved sheets. 



A major contribution of the NASA Structures Manual to the analysis of the skin and stringer 
structure is the determination of critical load that includes the load carrying capability of the post- 
buckled skin. From reference 8, the equation for critical buckling load (modified for integral stiffeners) 

is written: 


Per = (Fcr* J(A st H sk w e )+(F cr 

One can see the critical load in the column ( P C r ) is increased by the buckling stress of the 
skin multiplied by the area of skin not counted as effective skin acting with the stringer. However 
from an analysis procedure standpoint, it is easier to use the load carrying capability of .the buckled 
skin to reduce the ultimate stress in the column, rather than increase the column capability. 

Equation (28) is further modified by the observation that the skin does not carry the full buck- 
ling stress after failing. Bruhn 5 suggests that the maximum stress assumed in the buckled skin 
should be no greater than 


a cr = 0.3* E*tlR . 

In practice, the buckling stress value for the skin may be reduced 10 percent so that equation (28) 
becomes: 


p cr = (Fcr^J(A st +t sk w e H0.9(FcrJ(b sr w e ^ . (ref. 4) 


(29) 


The reader should note that the load carrying capacity of the buckled skin can be accounted for 
directly by raising the wide column buckling load or reducing the ultimate stress by reducing t e 
applied load, as has been discussed. Also, the buckled skin capacity for load can be retrieved 
indirectly by determining an area not included with the stringer or effective skin that will act to carry 
some load. This area, designated A e , will be used to increase the moment of inertia and total area 
calculations of the shell, and thereby, lower the stress level. Use of the ineffective area is covered 
under the section addressing stress calculations using simplified beam theory. The external tank 
(ET) stress method uses the P cr approach outlined above. 


2. Shear Buckling 

It is not often in launch vehicle design that shear buckling of a launch vehicle becomes a driv- 
ing load condition. In general, consideration of axial and bending loads far outweigh any consideration 
of shear effects. For completeness however, determination of the shear buckling capability of the 
skin will be presented here. 

The critical elastic shear buckling stress is given by the following equation: 


*cr 


(t\ 2 

12(1 -u 2 )UJ 


(ref. 5) 


If buckling occurs above the proportional limit, equation (31) must be employed: 


(30) 
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(ref. 5) 


(31) 


7 lAg 

cr 12(l-u 2 )Vf>) ’ 

where 7] s represents the plasticity correction factor. Correlation with test results indicates that an 
rj s = G/G yields best results. G s being the shear secant modulus, and G being the shear modulus. 
The coefficient k s is chosen from figure 14 using the hinged edge constraint (flat plate). 5 

Returning to the original reference by Gerard and Becker, a curve for shear buckling of a 
curved panel with simply supported edges can be obtained (fig. 15). 15 It is recommended that the 
largest of these buckling coefficients be used in determination of the shear buckling stress, keeping 
in mind that b is always the shorter of the panel dimensions. 
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3. Internal Pressure 

Internal pressure in a tank creates a biaxial tensile stress state in the skin and thereby 
increases its resistance to buckling. The addition of internal pressure increases the buckling strength 
of the curved sheet by the following interaction relationship: 6 

R 2 c +R p = 1. (32) 

where R c is the ratio of compressive buckling stress to critical compressive buckling stress, and R p 
is the ratio of applied internal pressure over the external pressure that would buckle the cylinder for 
which the curved panel is a section. Buckling due to radial pressure is found by use of figure 16. For 
internal pressure, R p is negative. The buckling equation is the same as for flat plate, substituting ky 

for k c . 3 * 5 
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Figure 16. Pressure buckling. 

A more sophisticated method of determining the increase in buckling strength due to the 
effects of internal pressure is found in reference 16. Assuming all edges are simply supported, the 
relation between critical meridional and hoop stresses is given by the following equation: 


^ + ^71 = 0.823-^t 

a 2 y b 1-u 2 



(33) 


Here m and n signify the number of half waves in the x and y directions, respectively. To find o y for a 
given g x , take m = 1, n = 1 if, 



where 

r 0.823E? 2 
( l-v 2 )a 2 

If g x is too large to satisfy the inequality, lake n = 1 and m to satisfy: 

C(2m 2 -2m+\+2 ,,2 /b 2 ) < G x < C(2m 2 +2m + 1 +2" 2 /h 2 ) . 
If G x is too small to satisfy the first inequality, take m = 1 and n to satisfy: 


C 


-n 2 («-l) 2 fr 


<G X <C 


\-n 2 {n+\) 


2 a M 
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( 34 ) 


Internal pressure also increases the shear allowable by the following relation: 5 

Rs+R p = 1 • 

This equation is employed in the same manner as equation (32). R s is the ratio of applied shear 
stress to the critical allowable shear stress for buckling. 

The alternative method presented by reference 16 maintains that the unit shear stress for 
buckling with all edges simply supported is given by equation (35): 



(a x and cr y are negative when tensile). 


III. EXAMPLE PROBLEM 


The launch vehicle shown in figure 17 will serve to illustrate the methods presented in this 
report. The hydrogen tank in particular will be examined. The diameter of the vehicle is 331 in. A 
Zee-shaped intermediate ring frame is shown in figure 18. The internal stringer geometry is shown 
in figure 19. 17 

The forces acting on the vehicle are due to ground winds acting against the vehicle prior to 
launch and the vehicle’s own weight. It is assumed for analysis purposes that no pressure is present 
in the tank at this time. At the hydrogen tank barrel section to be examined, the shear load is 
53,678 lb, the bending moment is 4.924x1 0 7 in-lb, and the axial compressive load is 1,571,825 lb. 
Recall that the skin is allowed to buckle at limit load. 

The skin and stringers are machined from aluminum 2219-T87 and the rings are extruded 
2024. A safety factor of 1 .4 will be used. 


A. Bifurcation Buckling 

A Fortran program written using MacTran — a Fortran development program for the 
Macintosh — was developed to facilitate the matrix method buckling analysis procedure. The code 
itself is found in appendix E along with a sample output. The code is not as autonomous as the 
worksheet developed for use with the other methods. Sectional properties for the shell configuration 
must be hard-coded into the Fortran program for execution. 
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Figure 17. Vehicle configuration. 


Figure 18. Intermediate ring frame. 
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Once again, recall that the bifurcation buckling methods may not be used in instances where 
significant redistribution of stress occurs in the subcritical load range — as occurs in this example 
problem. Determination of the bifurcation buckling allowable is presented for demonstration pur- 
poses. 


Using the same stiffened shell configuration shown previously, the stiffness matrix [Cy] is 
shown below: 


1,805,277.3 

503,943.5 

0.0 

-215,235.7 

0.0 

0.0 

503,943.5 

1,527,101.4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

511,578.9 

0.0 

0.0 

0.0 

-215,235.7 

0.0 

0.0 

210,811.6 

666.7 

0.0 

0.0 

0.0 

0.0 

666.7 

2,020.4 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

11,584.0 


N y is zero in this case since there is no pressure acting on the cylinder when the stated loads are 
applied. The minimum Nx load is found when m = 1 and n = 12. The classical buckling load is deter- 
mined to be 


N cl = 5,707 Ib/in. 

The wide column allowable is calculated as 

N wc = 1,137 Ib/in. 

The “knock-down” factor resulting from an (R/t) e value of 130.6 is 0.458. Therefore, 

N C r = 3,231 lb/in. 

In this particular case, the classical bifurcation load is of little use since there has been “consider- 
able redistribution” of stress (i.e., the skin has failed). 


B. General Instability 


The ring shown in figure 18 will be used as the intermediate stiffening ring in this example 
problem. The ring spacing is 40.08 in. Aluminum 2024-T42 extrusion has a compressive modulus of 
ll.OxlO 6 lb/in 2 and a compressive yield strength of 38,000 lb/in 2 . 11 


(3). 


Using Cf = 1/16,000, the critical equivalent bending moment can be solved for using equation 

(El)f = C f MD 2 /L , 


M eq = 2.7279xl0 8 in-lb . 


This equivalent moment is converted to a line load or stress for comparison to applied 
stresses in the structure. Comparison of the frame capability and the applied stress yields a margin 
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of safety of 0.428 using the Shanley method. See the section on margin of safety calculations for more 
details. 

Evaluating this frame using the Becker method the critical stress is 32,087 lb/in 2 . This allow- 
able stress determined by the Becker method equates to a margin of safety of 0.019. Evaluation of 
the (fi/Rt term reveals it to be less than 100; not in the range specified by Becker. As a test, a con- 
figuration was set up so that cP/Ri was greater than 100. For the test case, the two methods pro- 
duced identical results. 

The reader should take special note that the capability of the frames is compared to the 
applied stresses rather than to the equivalent moments themselves. The reason for this is that the 
stress in the shell is reduced by the load carrying capability of the buckled skin. This reduction in 
load is not reflected in the applied equivalent moment calculations, but is included in the Nx' value 
that is introduced as part of the ET stress method in appendix B, and in the final stress value pro- 
duced using the simplified beam theory. 


C. Stringers 


1. Local Elas tic Buckling 

The skin stringer panels under investigation are machined rather than formed or extruded. 
Divide the Tee stringer into three flange sections. Taking the conservitive definition of flange ele- 
ments, sections one and two represent the halves of the stringer cap. Section three is the web sec- 
tion, which actually extends to the centerline of the cap. 


in — 1 .250" ► 



The local buckling strength of flanges one and two, and that of flange three (the web) are 
determined in the following: 


F _ 0.43^ 2 1Q.8E6 
C| 12(1-0.33 2 ) 


(tif) 2 = l7 '- 452IM " 2 - 
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, _ 4.0;r 2 10.8E6 
Cj 12(1-0.33 2 ) 



= 367,584 Mn 2 . 


Obviously, the minimum buckling strength is above the material yield strength. Therefore, the 
local buckling limit is set equal to the yield strength of the material (51,000 lb/in 2 ). As the reader 
may have already discerned, calculation of the local buckling strength for this particular example 
problem is unnecessary. The stringers are internal to the tank structure and have no insulation or 
other bonded materials to be concerned with. The calculation is carried out for completeness and 
added confidence in the design. 


2. Cri pp lin g 

a. NASA Structures Manual . The following illustrates the procedures found in the NASA 
Structures Manual for determination of stringer crippling stress. 


Step (1) Maintaining an end-fixity coefficient of 1.5, L'/p equals 80.15. 

Step (2) Writing equations for the no edge free and one edge free curves from figure 7. 4 



1.387194 



'ficA _ 

Fey) ~ 


0.5693108 


OEF 



0.8127115 


The maximum cutoff value for each material is given by: 


feen ^tu 

Fey F cy 


or simply F tu • The procedures outlined make no mention of stresses in the inelastic region for 
stringer crippling. Obviously, if f ccn were to reach the cutoff limit it would be past the proportional 
limit stress and some correction factor would be in order. In this report, the cutoff stress will be 
limited to the yield stress. 



0.5693108 




51,000 

fO.575) 

10.8xl0 6 

10.1257 

i 


0.8127115 


1.45159 . 


The limit value is fccJFcy - 1.0; therefore, the crippling stress for the flange element is 51,000 
lb/in 2 . Notice from figure 6 that the flange width is not simply half the stringer cap width. 


Crippling strength for the web is determined using the no-edge-free equation as follows. 
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( fccn ) 1.387194 , o 1 

U^Lb"/ / 51,000 r 0.979 ^ a8g71793 ’ 

\V 10 . 8 xl 0 6 v 0.10 if 

Again the limit ratio of f ccn /F cy = 1.0 is employed and the crippling stress for the web 
becomes 51,000 lb/in 2 . 

The final effective width after iteration using the ET stress method is 4.28 in, so that: 

f fccn) 1.387194 i 2Q 

U^L"/ / 51,000 ( 2.09 I) 08071793 ' ’ 

\V 10.8xl0 6 10.1318// 

The resulting crippling stress for the effective skin is 51,000 lb/in 2 , and the weighted average 
of the three crippling stresses is also 51,000 lb/in 2 . Based on literature review and experience, it is 
recommended that the designer use the NASA Structures Manual method for determining stringer 
crippling first. The structures manual method is easier to program and more straight forward in its 
application. The analyst attempting to employ the Gerard method should review the papers authored 
by Gerard. The method of determining the number of cuts and flanges is somewhat confusing. 

b. Gerard . The shell geometry to be evaluated in this example is shown in figure 18. The 
Gerard method of stringer crippling analysis is obviously more applicable to the example than the 
Needham method. Dividing the Tee into angled sections, as called for by the Needham method, 
would be difficult without splitting the web. Therefore, the Gerard method is illustrated here. The 
Needham method will not be used in the example problem. 

Since the web and stringer cap are of different thicknesses, an equivalent thickness must be 
used. This 7^ is obtained from equation (17). 


w ~ Xbi ' 

To include the stringer land thickness with the rest of the effective skin, an average thickness 
is used. Figuring this average thickness in much the same way as 7 W ; t av = 0.1318 in. 

, 2(0.625)(0.125)+2(0.4895)(0.1)+2(2.14)(0.1318) _ n . 

(2(0.625)+2(0.4895)+2(2.14)) 

Then from equation (16), the stringer crippling stress with g equal to (7) (six flanges, one cut) is: 


F r , = 5 1,000(0.5346) 


7 (°. 1257)^T 13 18) ^do.gxto 6 ^ 1,000) 


0.85 


= 50,464 lb/in : 


P g is determined by interpolation from table 1. Table 2 shows the stringer crippling stress of a 
stiffened panel to be less than or equal to the proportional limit stress. Therefore, 


F cs = 50,464 lb/in 2 . 
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3. Column Failure 


Calculation of the column buckling stress involves the stringer properties as well as the 
effective skin width. Using equation (27) with the stress in the stringer equal to 31,503.6 lb/in 2 , the 
effective width is: 


w = Ut(E/F s ,) 0 - 5 =4.28 in. 


The properties of the stringer and the effective skin width are then combined to form a wide 
column. From appendix A, the stringer meets the requirements for a Tee section, and the effective 
width is simply the 4.28 in centered under the Tee web. The moment of inertia of the stringer and its 
effective skin is 0.131 in 4 . The end-fixity coefficient used in calculation of L' was 1.5, yielding an L' of 
32.7 in. Inspection of the L'/p value for this panel reveals that it lies in the Euler buckling region. The 
resulting critical wide column buckling stress is: 


F * 2 10.8xl0 6 
c (32.7/0.3997) 2 


15,926 lb/in 2 . 


Note that many values used in hand calculations are taken from the worksheet. Small differences 
may arise between those numbers and the ones shown here due to round-off errors not as prevalent 
in the worksheet. 


Had the L'/p value fallen in the transition region, the Johnson-Euler equation would have 
been applied as shown below: 


F c 


= 51,000- 


(51, OOP) 2 f 32.7 ) 2 

4;r 2 (10.8xl0 6 ) v0.3997/ 


10,169 lb/in 2 . 


It is useful here to review the transition from Johnson-Euler to Euler equations. This transi- 
tion occurs at the intersection point of the two curves when plotted as functions of L'/p. This critical 
L'/p value is found by setting the two equations equal to one another. If the L'/p value is larger than 
the critical value, use Euler. If the L'/p value is smaller than the critical value, use Johnson-Euler. 


( L 7 P)critica]=*y ~jr ■ 
In this case, the critical value is L'/p = 64.65 in. 


D. Skin 


1. Compression Buckling 

For compression buckling of the skin, equation (26) is applied with b equal to the stringer 
spacing and t the skin thickness. The buckling coefficient value most used by Bruhn and others is the 
conservative 4.0. The compression modulus for A1 2219 is 10.8xl0 6 lb/in 2 , and u is 0.33. 11 
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12(l-u J )UJ ' 

Using 4.0 for k c yields a critical stress of 5,395 lb/in 2 , a very low number indeed. Buckling 
stresses are typically well below the material limits (F lu 2219-T87 = 63,000 lb/in 2 ). If a less conserva- 
tive approach is taken and k c is chosen from figure 12 (k c = 5.80), the critical stress increases to 
7,823 lb/in 2 — a considerable increase. 

Using data for curved simply supported sheets, and reading k c from figure 13 with Z = 5.31 1, 
the buckling coefficient is 5.09. The buckling stress is then computed to be 6,861 lb/in 2 . Note here 
that the buckling coefficient, and therefore, the buckling stress for the curved panel, is lower than that 
of the flat sheet. The reason for this apparent contradiction is that the curved panel data is for simply 
supported edges, while that for the flat sheet is for an edge condition between simply supported and 
clamped. It is recommended that the maximum of the flat sheet and curved panel buckling coefficients 
be used. Therefore, the critical compression buckling stress for the skin between the stringers is 
7,823 lb/in 2 . 


Shear buckling is handled much the same way as compression buckling. The buckling coeffi- 
cients are read from figures 14 (flat sheet) and 15 (curved panels). The buckling coefficient for flat 
sheet is 5.8, while for curved panels it is read as 6.05 (simply supported curved panels). The result- 
ing shear buckling stress is 8,160 lb/in 2 . Again, the maximum is used and the critical shear buckling 
stress is reported as 8,160 lb/in 2 . 


E. Example Summary 

Presented, thus far, in this report have been the basic methods of stringer stiffened panel 
design as presented by Dr. Bruhn, Mr. Almroth, and others. The methods used have been gleaned 
from many portions of Bruhn’s books, books and papers by Almroth, and various other government 
and journal publications. It should be noted that the methods presented by Bruhn were devised 
before the advent of modern computational devices. Therefore, most of the methods rely on design 
curves to lessen the computational intensity of the problem. However, curve fitting of the appropriate 
design curves can make the methods acceptable to modern programming techniques. All necessary 
design curves have been curve fit and programmed as macro routines in the Excel programs 
contained in the appendices. Despite their age, many of the methods compiled by Bruhn and 
presented in this section are still used extensively in the aerospace structural design field. An 
engineer must understand and know how to apply these methods before exploring the more recent 
works. A summary of analysis results is shown in table 3. 
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Table 3. Summary of critical stresses. 


Failure Mode 


Critical Stress (lb/in 2 ) 

General Instability 

Shanley 

45,000 


Becker 

32,087 

Local Crippling 


51,000 

Column Failure 


15,926 

Skin Compression 


7,823 

Skin Shear 


8,160 


The following section of this report details the determination of applied stress and margins of 
safety. Obviously, the applied stress cannot be determined independently of certain failure modes 
already covered in detail. Discussion of applied stresses and margins of safety is placed in a 
separate section for organizational purposes only— to separate failure analysis techniques from 
applied stress calculations. 


IV. APPLIED STRESSES AND MARGINS OF SAFETY 


This section of the report provides a summary and explanation of loads and stresses applied 
to the example configuration. Determination of the compressive stress level at limit load is fairly 
simple — all the skin is effective and moment of inertia and area calculations are quite simple since 
they involve the stringers and the entire skin rather than portions that are “effective” and portions 
that are “ineffective.” Equations (36) and (37) are used to determine the maximum compressive 
stress level. 

= , (36) 


/T -£ 

a ox~ A 


(37) 


Determination of compressive stress at ultimate load can be considerably more involved. If 
the skin has not failed at ultimate load, the stress calculation procedure is identical to that for limit 
stresses. If the skin has failed however, the process becomes a bit more cumbersome. Ultimate 
stresses in the shell can be calculated through several methods with varying degrees of accuracy. 
The most accurate calculation, of course, requires the most rigorous analysis of the shell configura- 
tion. If the skin buckles after limit load, the effective skin provides some stress relief to the stringer 
columns. This stress relief depends upon the effective width of skin acting with the stringer. The 
effective width is dependent on the stress. One can easily visualize the iterative process necessary 
for determining stress and effective width. Shear stress calculations are performed at the limit load 
level and the procedure is quite well known. Calculation of the shear stress is shown in detail in 
appendix B. 

Figure 18 gives details of the configuration to be evaluated. 17 The worst case mechanical 
loads (limit loads) on the shell are shown in the following and are typical of launch vehicle ground 
wind induced loads. 18 
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Shear = 53,678 lb 

Moment = 4.924xl0 7 in-lb 
Axial = 1,571,8251b. 

This particular load case occurs during the prelaunch phase, and the vehicle at this point is con- 
sidered as being fully fueled and unpressurized. 


A. Simplified Beam Theory 

The most accurate determination of stress involves calculation of an initial stress estimate 
using the entire skin as effective. Since the bending moment acts in compression on one side of the 
neutral axis and in tension on the other, the neutral axis of the shell will be shifted toward the tensile 
moment side of the shell. The neutral axis shift is a direct result of increased effective width at the 
lower compressive stress. 

Figure 20 illustrates the iterative shift in neutral axis caused by the application of bending 
moment. 5 



Figure 20. Illustration of neutral axis. 


The procedure begins with an estimate of stress using the moment of inertia and area includ- 
ing the entire skin. At this point the stress in each bay must be calculated individually along with an 
accompanying effective width. When this is done, a new area, moment of inertia, and neutral axis is 
computed using only the stringer, its effective skin, and A e as defined below. Including this area (A e ) 
can be viewed as having the same effect as the reduction in load presented by the NASA Structures 
Manual. Stress in each bay is then computed with the new cross-sectional properties. This proce- 
dure is repeated until the neutral axis location converges. 


K = (K r l Mw e . ,-w )) 


CR, 


(ref. 5) 


(38) 
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The symbols w e i and w ei+ i refer to the effective widths on either side of the stringer location where 
Ae is desired. Since the effective width changes very little between adjacent stringers on a large 
diameter cylinder, equation (38) can be rewritten as equation (39): 


A (U \\ a<:R i 

A e=(b s r”e)[-0^ • 

The stress in each bay may be calculated on the stringer location by: 

My: j P 

a b=— and ^ax=j: ■ 

The y, in this equation refers to the stringer distance from the centroid. 


(39) 


B. ET Stress Report 

This method is referenced from the ET stress report produced for NASA by the Martin 
Marietta Corporation. The method in the stress report lacks a great deal of referencing, but has been 
used successfully in the analysis of the Space Transportation System ET. The main advantage of 
this method is that it is much simpler, requiring iterations only with the maximum stress in the shell. 
That is, stresses at each station about the circumference of the shell do not have to be calculated. 
This greatly reduces the computational intensity of the process. Also, this method allows inclusion of 
hoop stresses in determining the longitudinal stress. This feature could be included in the simplified 
beam theory method with a bit of derivation. The disadvantage of the method is its conservatism. 
Choosing the proper method depends on the analyst’s expectations of fidelity in the analysis, confi- 
dence in the given load set, and ultimately, the cost of failure. 

Using this method, a more detailed stress breakdown can be obtained for the skin, stringer, 
and stringer land. The equations for normal stresses in the skin, land, and stringer are presented 
along with figure 21 for explanation. 4 



Figure 21. Stress distribution in shell. 
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* sm 

fld = T~~ 


N x +vN y 


N x +vN y 


(%?)]• 

m- 


f,,= 1 L H- vN >] ■ 

l sm L 


The line loads (N y and N x ), in pounds per inch, are determined based on stresses developed 
in a thin cylindrical shell by internal pressure and axial and bending loads, respectively: 

,r _ n si _ Axial j_ Moment 

N y = pR and N,__±— p- , 

and t sm is the smeared thickness of the stringer, land, and skin thickness so that, 

t 

l sm U 

u st 

Stress in the stringer is now determined based on the applied line loads and the shell 
geometry. Now, providing the skin buckles at F cr , the longitudinal load capability of the skin panel is 
0.9 *F cr *t skj or 90 percent of its buckling load. This capacity of the buckled skin to carry load reduces 
the compressive load that the stringer must support, just as A e adds to the load carrying capacity of 
the shell. The effective width acting with a stringer is calculated by equation (27). This varies 
somewhat from the ET stress report method for determining effective width, but should give com- 
parable results. Once the effective width is known, the load carried by the ineffective skin is deter- 
mined and subtracted from the total load applied. Line loads are converted to forces in each panel by 
multiplying by a characteristic length — either stringer spacing or effective width depending on 
whether or not the skin has failed. The new load that must be carried by the stringer is now known 
and the process repeats itself until convergence is obtained on the effective width or the stress. 


C. Margins of Safety 

Table 4 is a summary of stresses produced using the ET stress method and simplified beam 
theory. As shown in the stress table, the shear stress is very low compared to the compressive 
stresses resulting from bending and axial load. Shear stress is often ignored for preliminary analy- 
ses. 

Table 4. Applied stresses. 



Limit (lb/in 2 ) 

Ultimate (lb/in 2 ) 

Shear 

Moment 

Axial 

Total Compressive 

819.36 

3.770.4 

9.959.5 
13,729.9 

31.503.6 (ET) 

24.645.7 (Beam) 
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The margin of safety (MS) is a numerical evaluation of a structure’s load carrying capacity 
compared to the applied load. In general, when there is only one type of loading, the MS calculation 
takes the following form: 5 


_ Allowable Stress or Load _ j q 
A pplied Stress or Load 


(40) 


Evaluating the general instability equations first, the critical stress determined using the 
Shanley method is 45,000 lb/in 2 . However, 45,000 is greater than the yield stress for the ring of 
38,000 lb/in 2 . Therefore, the capability of the frame is limited to the yield stress. The applied ultimate 
stress is 31,503.6 lb/in 2 . The following MS calculation results: 


ms = - 3T1M _1 - 0 = 0 - 2062 ' 

For the stringer and its effective skin, the wide column buckling allowable is lower than the stringer 
crippling value. Therefore, the wide column value is used to determine the minimum MS. The applied 
compressive stress is 31,503.6 lb/in 2 , and the critical buckling stress is 15,917 lb/in 2 . The resulting 
margin of safety is -0.495. (The shear load is ignored for the wide column margin of safety because of 
its relative insignificance.) Obviously the column is inadequate for the applied loads and must be 
redesigned. 

The skin is under combined shear and compressive load. The margin of safety under combined 
shear and compression loading from reference 5 is, 


MS = 


R c + s/ Rc+4*R 


’ — 1 


(41) 


where R c is the applied compressive stress divided by the critical buckling or allowable stress. R s is 
the applied shear stress divided by the allowable shear stress in the skin. The MS for the skin is 
determined using limit loads. The resulting MS for the skin is -0.4316, with R c = 1.7551 and 
R s = 0.0874. The shear stress contribution to MS is practically negligible. This is typical in launch 
vehicle design. As with the column buckling MS, the skin is inadequate for the loads applied and 
must be redesigned. 


V. CONCLUSIONS 


In this report, many of the most popular methods for determining buckling capability in a ring 
and stringer-stiffened cylindrical shell have been presented. Methods for determining the skin buck- 
ling load, the stringer failure allowable, and the wide column allowable have been presented. Two 
methods for determining the necessary ring geometry to preclude general instability have also been 
presented. Where applicable, the conservative approach taken by most designers has been pointed 
out, along with methods for reducing unnecessary conservatism. 

Microsoft Excel spreadsheets have been developed in conjunction with this report to facilitate 
the use of the methods presented. All design curves necessary for calculation of critical buckling 
allowables have been curve fit and included in the worksheet as macros which act as subroutines for 
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calculation of buckling coefficients from the curves. Calculation sheets done by hand have also been 
included in the appendices. The most common problem facing the analyst in this area is the proper 
designation of sectional properties, i.e., when to include effective widths and when not to. 
Considerable effort has been made to track down any discrepancies of various reports on the use of 
sectional properties. Where clear instructions were not given by the author, the most logical option 
regarding the use of sectional properties was chosen. 

Although limited in its scope of application, a method for determination of the bifurcation 
(classical) buckling load has been included. This method cannot be used in cases where considerable 
redistribution of stress occurs in the subcritical load range — as happens with the example problem 
presented. A Fortran program has been written to facilitate determination of the buckling allowable 
using the methods devised by Almroth. 

The report also addresses determination of the applied stress and MS. Neither the applied 
stress nor the buckling capability of the shell can be determined independent of one another. 
Recognize that the MS is the measure of structural margin most often used in conjunction with the 
safety factor. 

This report has been designed to serve as a reference for the analyst in need of determining 
the buckling capability of stiffened cylindrical shells. The methods are easily modified for use with 
stiffeners other than the integral Tee’s used for the example analysis. However alterations to the 
programs would be necessary. Critical stress calculations for the various methods are given in the 
following chart: 


Summary of Critical Stresses 


Failure Mode 


Critical Stress (lb/in 2 ) 

General Instability 

Shanley 

38,000 


Becker 

30,930 

Local Buckling 


51,000 

Crippling (NSN) 


51,000 

Crippling (Gerard) 


50,464 

Column Failure (Euler) 


15,917 

Skin Compression 


7,823 

Skin Shear 


8,160 


The methods presented in this report are not meant to exclude or replace the use of analysis 
codes such as BOSOR or PANDA. Neither have all possible failure modes for tank design been 
addressed. The report has addressed the major failure modes associated with stiffened circular 
cylinders and provides the methods necessary for assessing vehicle design. Once the major 
geometric properties of a design have been established, more rigorous analysis using further refined 
hand techniques and computer algorithms would be in order. 
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Also not covered in this report are the optimization techniques. Considerable work has been 
done in the optimization of stiffened shells area. Most of the methods are mathematically derived, 
rather than being based on empirical data as most of Bruhn’s methods are. It is suggested that the 
engineer first become familiar with the methods presented in this report, and then investigate the 
various optimization techniques. An optimized design will likely need to be checked against the 
methods presented in this paper to ensure the capability of the design. 
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APPENDIX A 


Tee Stringer Criteria 
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APPENDIX B 


Hand Calculations 


General Instability 


Shanlev 


Becker 


m f =c,MDi 


[Bruhn, "Aircraft" equation (9.7)] 

Cf= 1/16,000 


Solving for M => M CR = 


m f L 

C f D 2 


FRAME = 2024 A1 

Ef= 11.0E6 in-lb 
I f = 4.2369 in 4 


Mcr - 


( 1 1 .0E6)(4.2369)(40.08) 
(1/1 6,000)(33 1 ) 2 


M cr = 2.7279x 10 8 in-lb 


(o. o 


i 


Fcr = CE c Qf, 

[NACA TN 3786: pg. 33] 

( Ps p f ) m W 25 


Q„ = 


R 


C = 3.25 Frame not attached to skin 
C = 4,80 Frame attached to skin. 

p s = radius of gyration of stringer and effective skin 


W' 


w. 


„ = °.5 5fiR7-(l^) 


10.832 


determined by: = 0.5 SQRT 

Fqr = Buckling stress of column 


m 


H’ =4.41 in 

e & 


F st = Stress in stringer 
A effa = 0.810 in 2 =0.171 in 4 P f = radius of gyration of frame and effective skin. 

w e f = d 


Pst = 


= 0.4595 in 


Qb = [(0.4595)(1.3789)Qq0.832)(40.08)]- 1/4 = 0 000940 

F CRm = (3.25)(1 1.0£6)(0.00094) = 33,608.12 lb/in 2 . 

*Note that moment of inertia of frame is computed without the large stand-off caused by mounting on 
top of the intergral stringers. 
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b = 1. 104 = 1.0415 in 


_ 4.0 k^_ (10.8E6) f Q.j_^ _ 367,584 lb/in 3 
12(1-0.33 2 ) v 1.0415; 


Minimum Fqr > Fcy Fqr - Fcy 
Local Buckling Limit = 51,000 lb/in 2 
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1.104 



Stringer Cripp lin g 


♦ NASA Structures Manual 



U)e > 



1.387194 


NO EDGE FREE 


FcyY^bAV ^ 793 




0.5693108 


f ONE EDGE FREE 


Fcyl^fbA) 08127115 


4 M 

Ec\ Vn)! 


Curve fits from NASA 
Structures Manual Cl. 3. 1-13 


Corresponding to those 
in ET stress report. 


(44 

0.5693108 

144 

/[ 5 1,000 ] m (0.575 08127115 


\[ 10.8E6J 10.125// 


= 1.4516 


in 2 


f cc = (1.4516X51,000) = 74,030.83 > 51,000 /. 4 =51,000 lb/in 


(44 

144 

( Xc 3 ^ _ 

144 


1.387194 


FcyV a 


(0.979 u 
l 0.1 )J 

1.387194 


0.8071793 


= 1 .9 1 > 1.0 fcc 2 = 51 .000 lb/in 


H/2 


O' 


2.09 1 
0.1318J 


0.8071793 


= 1.29 >1.0 


^av 


4.28(0. 126)+1. 25(0.02) 
4.28 


4 = 51,000 


^=0.1318 
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Stringer Crippling 


'LbMc, 


[NASA Stress Manual, Section Cl, page 11, equation (1)] 


F cc = 5 1,000 { (0.575)(0. 125)(2)+(0.979)(0.1)+2*(2. 14X0. 13 1 8) }+ 

{ (0.575X0. 1 25)(2)+(0.979)(0. l)+(2. 14X0. 1 3 1 8X20) } 


Stringer Cripp ling 


Gerard : 


w e = 4.28 , bf= 1.25 , h w = 0.979 
Using an average thickness to evaluate the effective skin with the land: 

, 1.25(0.02)44.28(0.126) A ,„ 10 

* av~ 4 28 — U.lJlo 

_ _Zb i t i 2(0.625)(0. 125)+0.979(0. l)+4.28(0. 131) 
* w ~ Y,b, " 2(0.625)+0.979+4.28 


7 W = 0.1257 in 

4 = 1.25 (0.1 25)+0.979(0. 1)+4.28(0.13 18) = 0.8183 in 2 


=> Interpolate 


_ 0,1257 
t sk "0.1318 


0.954 


1.16-0.732 0.562-0.505 

0.954-0.732 " /J-0.505 

1.9279 = 0-ffM £-0.505 = 0.0296 
0 = 0.5346 



(7) 


wX^qv) 

4 



1/210.85 


F„=0.99/v, . 

F w = 50,463.89 lb/in 2 = 50,464 lb/in 2 
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Column Buckling 


Determine Johnson-Euler Intersection 



[2 E c ] 

f L'\ — n 

l P / INTERSECTION 

.Fes. 


( (2)(10.8 E 6)\ 1/2 

\ 51,000 / 


= 64.65 in 


_(I\ m 

P COLUMN ~ \ Af 


I and A = Sum of stringer and effective width. 


w e = 4.28 in 



■^COLUMN = 0.254+ 1 .25(0.02)+frfcWe 
^column = 0-82 in 2 


Centroid = 


SM 

Za, 


Centroid = [( 1 .25)(0. 1 25)( 1 . 1 875)+(0. l)(0.979)(0.6355)+( 1 .25)(0.02)(0. 1 36) 


+ w e ^(f^/2)]+4coLUMN = = 0.35 in 

/column^ ^ (1-25X0.125) Vl.25(0.125)(l. 1875— 0.33) 2 +i( 0 . 1 )(0.97 9 ) 3 +( 0 .l)(0.979)(0.6355^.33) 2 

+ i ( 1 .25)(0.02) 3 +( 1 .25)(0.02)(0. 1 36-0.33) 2 + ± w e t] k +w e t s ^k/ 2 -0.33) 2 

^column = 0.131 in 4 


P = 



0.131 

0.82 


= 0.3997 in 


Lj 40.08/vT5 
p ~ 0.3997 


81.87 
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Column Buckling 


Using an end-fixity coefficient of 1.5 and a ring spacing of 40.08 in, the effective length is 32.72 
in. When divided by p, this length becomes the slenderness ratio. 


^ = 81.83 > 



INTERSECT 


Column is in the Euler regime. 


Critical buckling load for the column is: 


F _ ;r 2 E = ;r 2 (10.8E6) 
CR (81.83) 2 


15,918.3 lb/in 2 . 
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Skin Buckling 


Compression Buckling of Elat Sheet 


Fcr - 


12(l-v 2 ) 



[Bruhn, “Aircraft,” equation (C5.1)] 


From figure C5.6 with j = = 85.968 

k c = 5.802. Could use 4.0 for conservatism. 


F cr = 5.802 ;r2(10 ' 8£:6 2 ) f-M^) 2 = 7,825.7 lb/in 2 
12(1-0.33 2 ) v 10.8327 


Compression Buckling of Curved Panel 


k c n 1 E c ( t \2 
CR 12(l-v 2 ) V6/ 


From Bruhn, “Aircraft ” figure C9.1 Z = ^ (l-v 2 ) m 


Z = 


10.832" 


(165.5K0.126) 


From figure C9. 1 


(1-0.33 2 ) i/2 = 5.311 with R/Z = 1,313.49 
k c = 5.09 


r- 5.09 nj 10.8E6 
CR 12(1— v 2 ) 


(■ few ) 2 = 6,86 °* 5 lb/in2 


Taking the larger — Fcr = 7,825.7 lb/in 2 
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Skin Buckling 


Shear Buck ling of Flat Panels 


From figure C5. 1 1 k s for 


r - k ^ 2R 

Z CR~ 


' r. / \ 2 

(^J [Bruhn, “Aircraft,”” equation (C5)] 


12(1— v z ) 


a _ 40,08 _ t 7 _ o 
b~ 10.832 -i-'-™ 

k s = 5.8 (hinged edges) 

1 2(1-0.33 2 ) tali -^2.S9Mm- 
For simply supported curved panels . 

k s = 6.49 

r CR = 8,753 lb/in 2 . 
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A pplied Compressive Stress 


• Stress at Limit Load: S.F. = 1.0 

Limit loads are: moment = 4.924x 1 0 7 in-lb 

axial = 1,571,825 lb 

_ _ mz - _ P. 

°b-— <Ta ~ A ’ 

Ob = bending stress 
o\ = axial stress 
m = moment 

z = distance from centroid 
I = moment of inertia 
P = axial load 
A = area . 

*At limit load, all skin is effective. Therefore, 

I = (number of stringers) ; let t sk = t to include land 

/ = 96*0.0277+0.254 XZ 2 +;r(165.5) 3 7 

_ (0.02)(1.25)+(0.126)(10.832) jr 

* ~ 10.832 

7 = 2,161,079.1 in 4 

Area = InR't+lAsr = 2 tc( 165 . 5 )( 0 . 1283 )+( 96 )( 0 . 254 ) 

A = 157.79 in 2 

(4.924x10 7 )(165.5) 

Maximum bending stress = <T ^ = 2,161,079 

o b =3,770.9 lb/in 2 

^max 

. 1,571 ,825 o q'i 

Axial stress = O A = j 57 79 ” = a a ~ 9,960.93 


Total compressive stress = 0i, max +0/i 
o T = 13,731.83 lb/in 2 
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A pplied Compressive Stress 


• Stress at Ultimate Load: S.F. = 1.4 

Ultimate loads are: moment = (4.924xl0 7 )*1.4 in-lb 

axial = (1,57 1,825)* 1.4 lb 


Check Skin Buckling: 


cr. 


= _k£%_f 1 y 

12(1 -v 2 )\b) 


[ref. Bruhn, “Aircraft,” (C5.1)] 


k c = max from curves C5.6 or C9 . 1 


k c = 5.8 


_ 5.8/r 2 (10.8E6) ( 0 126 j 2 _ 7 ,822.89 lb/in 2 
CR * 12(1-0.33 2 ) 1 10.832; 


If skin does not buckle, / and A remain the same as limit. 

(4.924xl0 6 )(1.4)(165.5) 


o h = 


2,161,079 


= 5,279.26 lb/in 2 


13.945.3 


a T = 19,224.56 lb/in 2 

13,945.3 > 7,822.89: • Skin buckles prior to ultimate load causing redistribution of stress. 

* Note that <T ALjmit also > 7,822. This indicates a negative margin for skin buckling. The skin must be 
redesigned to prevent buckling at limit load. 
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A pplied Compressive Stress 


• The skin fails prior to ultimate load causing a redistribution of stress. The stresses at ultimate load 
must now be determined based on the effective widths acting with each stringer. 

* Simplified Beam Theory 

Step 1 : An initial estimate of stress in each stringer bay is made based on the entire skin as 
effective. 

7 = 2,161,079 in 4 

A = 157.79 in 2 


m = (4.924x1 0 7 in-lb)*(1.4) 


P 


Shell Geometry 



z\ =0 

Z 2 = sin(3.75) 165.5 = 10.8242 in 
z 3 = sin(7.5) 165.5 = 21.60 in 
z 96 = -10.8242 
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A pplied Compressive Stress 


* Note: A e 

b'-- 

* Note: A e i 


1,571,825*1.4 no ,, noi ,, 2 
a A = = 13,946.09 lb/in 2 

„ _ n tr (4.924x 10 7 )(1.4)(10.8242) ^Slh/in? 
^6,-0 , (7 b 2 ~ 2 161 079 — 345.28 Ib/in 


° b . 


(4.924xl0 7 )(1.4)(2l.6) _ 2 

2,161,079 - 689.02 lb/in 




(4.924xl0 7 )(1.4)(-10.8242) _ 2 

i" 2,161,079 " 345.28 lb/in 


<t , = 1 3,946 lb/in 2 ; w ti = 1 .7 (t s +t f )l2 
[ref. Bruhn, “Aircraft,” equation C7.16] 

w ei = 1.7(0. 1360) = 5.9608 in 

A e = ( 1 0.832-5.9608)(0. 1 26)[ 7 | ; ] 


A, =0.3443 in 2 


c 2 = 13,946+345.28 = 14,291.28 lb/in 2 


w, =5.8884 in , A e =0.3410 in 2 


c 3 = 13,946+689.02 = 14,635.0 
w e =5.8188 in , A e =0.3376 in 2 




ct 96 = 13,946-345.28 = 13,600.72 
w e = 6.036 in , A e =0.3476 in 2 

^96 € 96 

= b't(F c /<j,) (ref. Bruhn, “Aircraft,” equation (2), page A20] 

= bsT-w, 
s not w e *t sk \ 


58 


Applied Compressive Stress 


Step 2: Determine new centroid using A e , w e *t aVi , and Asj. 

(0.02)(1. 25)+(0.126)(vtv) 
t av . is new t ^ 

e i 

New centroid = -11.58 in 

New moment of inertia = 1,859,548 in 4 



Step 3: Recalculate stresses, A e and w e based on new sectional properties. 

(4.924xl0 7 )(1.4)(0+l 1.58) , (1,571,825)04) 
1,859,548 133.6 

<7, =429.28+16,471.2 = 16,900.5 lb/in 2 


(4.924xl0 7 )(1.4)(10.8242+11.58) . 

= 1,859,548 + 1M/LZ 

(T 2 = 17,301.7 lb/in 2 

(4.924xl0 7 )(1.4)(- 10.8242+ 11.58) . 1jr 0 
1,859,548 + lo,4/ t.z 


A pplied Compressive Stress 

Step 4: With new stresses, A e and w e , recalculate neutral axis location, moment of inertia, and area. 
Step 5: Repeat until convergence on neutral axis location is obtained. 

Final Results: 

Moment of inertia = 1,713,065.85 in 4 

Area = 124.185 in 2 

Neutral axis or centroid = -1 1.654 in 

<W = ct 25 = 24,645 lb/in 2 
w e25 = 4.84 in. 
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A pplied Compressive Stress 


External Tank Stress Method 

Step 1: Determine initial line loads. 

Ar 1,571,825*1.4 (4.924 x10 7 )(1.4) 

2 kR + *r 2 

= 2,917.3 lb/in 

N y = (P ress " re )( R ) . p ressnre = o,N y = 0 
Step 2: Determine initial stress based on line load and smeared thickness. 

, 2,917.3 

Jsk t 

* sm 

. _ - 0 - 254 + 1 . 25 ( 0 . 02 )+^/^ 


t sm - 0.1518 in 


j ! * = ^ W =19 ' 223 ' 49lh/in2 

Step 3: Calculate effective width. 

w e = 5.48 in . 
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Applied Compressive Stress 


Step 4: Determine new loads in the column. 

The total equivalent axial load carried by the panel is. 

Peg = 2,917.3*10.832 = 31,600.2 lb 
The load carried by the buckled skin is, 

P sk = 0.9 F cr t sk (10.832-wv) ; w e = 5.48 
P sk = 4,764.9 lb . 

Load which must be supported by the column is 

31,600-4,764.9 = 26,835 lb. 

The resulting line load, t m , and stress are: 

N x '= 26 ,f 35 =4,896.9 lb/in 
W € 

> 0.254+1. 25(0.02)+r Jt w e 

Lm “ Wg 

t'sm = 0- 1^69 in 

jS, ' = lw =27 ' 68a36IWi " 2 - 

Step 5: Continue iteration until convergence on stress is reached. 

The final results are: 

w e = 4.28 in 
hrn = 01912 in 

N; = 6,020.26 Ib/in 

f sk = 31,488 lb/in 2 = 31,503 lb/in 2 from spreadsheet. 
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Spreadsheet Output 


Skin-Stringer Spreadsheet 
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Skin-Stringer Spreadsheet 
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Fty = 51 OOOpsi Ec = 10.8E+6psi Fty= 38000 psi Ec= 11.0E+6psi 

Fey = 51 OOOpsi Fey = 38000 psi 

density = .10 lb/in^3 density = .100 Ib/irV^ 

Poisson= 0.33 _____ - 










Control Macro 


sesame (o) 

=OPEN("MOI“) 

=HIDE() 

=OPEN('CLOSER') 

=HIDE() 

=OPEN(’BruhnFig') 

=HIDE() 

=OPEN('l interpolate', 1 ) 
=HIDE() 

=OPEN(’Stress',1) 

=HIDE() 

=OPEN('Skin-Stringer') 
=ACTIVATE("Control Macro 0 ) 
=HIDE() 

=RETURN() 
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MOI 


Moil 

=ARGUMENT("Rad") 

=ARGUMENT("Nost") 

=ARGUMENT("Ast“) 

=ARGUMENT('lsr) 

=ARGUMENTnsk") 

=ARGUMENT('bf") 

=ARGUMENT(tl") 

bst=2*PI0‘Rad/Nost 

DA=2*PI0/Nost 

Ald= (tl+tsk)*bf 

tav= (tsk*(bst-bf)+Ald)/bst 

ltotal=0 

moist=0 

moisk= PI0*Rad A 3*tav 
=FOR('CountM,Nost) 
d=Rad*SIN((Count-1 )*DA) 
=SET.NAME(*moisf, moist + Ast*d A 2) 
=NEXTQ 


=SET.NAME(’ltotar,moisk+moist+Nost*lst) 

=RETURN(ltotal) 

Moi2 

=ARGUMENT('Rad“) 

=ARGUMENT('Nost“) 

= ARGU M ENT (' Aeff ") 

=ARGUMENT(‘leff“) 

=ARGUMENT(tsk“) 

=ARGUMENT( , bf") 

=ARGUMENT(’tr) 

DA=2*PI()/Nost 


ltotal=0 

moist=0 


=FOR(“Count",1 ,Nost) 
d=Rad*SIN((Count-1 )*DA) 
=SET.NAME(“moist“,moist + Aeff*d A 2+leff) 
=NEXT() 

=SET.NAME("ltotal", moist) 
=RETURN(ltotal) 
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CLOSER 


Closer (c) 

=SAVE() 

=CLOSE0 

=UNHIDE('BruhnFig") 

=SAVE() 

=CLOSE0 

=UNHIDE( , lnterpolate'') 

=SAVE() 

=CLOSE0 

=UNHIDE("Control Macro") 
=SAVE() 

=CLOSE0 

=UNHIDE('MOI') 

=SAVE() 

=CLOSE0 

=UNHIDE('Stress‘) 

=SAVE() 

=CLOSE0 

=UNH I DE("Closer") 

=SAVE() 

=CLOSE0 

=RETURN() 
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BruhnFig 


FigC5.6 

= ARGUMENT ("BOT") 

= ARGUMENT ("Zappa") 

=2.9384+0.088965*BOT-0. 001 3486*BOT A 2+0. 0000 12063*BOT A 3-0. 0000000531 53 *BOT A 4+0.000000000089722*BOT A 5 TorsionaHy Weak 

=4.1 134+0.074296*BOT-0.001 1 77 * BOT A 2 +0.00001 1 18*BOT A 34.000000052865*BOTA4+0.000000000095262*BOT A 5 TorsionaUy Strong 

=IF(AND(AND(BOT>»l5,BOT<=200),Zappa=1),RETURN(A4)) 

=lF{AND(AND<BOT>=10 1 BOT<=145),Zappa=2),RETURN(A5)) 

=IF(BOT<0,RETURN("bA Out of Range.")) 

=IF(AND<BOT<15.Zappa=1) f RETURN(4)) 

=IF(AND<BOT>200 r Zappa=1 ), RETURN<6.96)) 

=IF(AND(BOT<10 f Zappa=2),RETURN("bA Out of Range")) 

=IF(AND(BOT>1457appa=2).RETURN(6.96)) 


=RETURNfMacro Err") 


FigC5.11 

= ARGUMENT f AOB") 

=121.37-3t4.43*AOB+-350.03*AOB A 2-195.9*AOB A 3+*54.702*AOB A 4-6.0732*AOB A 5 
=-1 90.31 +368. 82* AOB *261 .64*AOB A 2+91 .472* AOB A 3*1 5.0* AOBM+ 1 .0804*AOB A 5 
=6.5950-0 -20958* AOB 


^<aAx22 

2.2<artx3.68 

3.68<a/t><5.0 


=IF(AND<AOB>=1 ,AOB<2.2),RETURN(A20)) 
=IF(AND<AOB>=22,AOB<3.68).RETURN(A21)) 
=IF(AND(AOB>*3.68,AOB<=5) t RETURN(A22)) 
=IF(AOB<1 ,RETURN(’a/b Out of Range")) 
=IF(AOB>5.RETURN("a/b Out of Range")) 

= RETURN ("Macro ErO 


RgC9.1 

=ARGUMENT(T') 

= ARGUMENT("ROT") 

=3.8337+0.25748*2-0.0015272*2*2+4.0000048691 *Z A 3 rA=3000 

=3.8337+0^5748*2-0.001 5272*Z*2 +0.0000048691 *Z A 3 rA=2000 

=3.4625+0.3351 *Z-0.0061 366*Z*2+0. 000084 875* Z A 3 r *= 1 000 

*5.5977-O.087272*Z+O.0l 5305*Z A 2 4.000 1 4876*Z A 3 M=70 ° 

=4.121+0.078764*Z+0.02733rZ A 2-0.001383*Z A 3+0.00003126rZ A 4-0.00000024929*Z A 5 rA=500 

»5.0601-1.7512*2+1.0254*Z A 2-0.24468*Z A 3+0.027988*Z A 4-0 0012329*Z A 5 Z>M AND Z<=7.0 

-A Z<=1 4 


=IF(AND(Z>= 1 ,Z<=1 .4),RETURN(A40)) 

=1F(AND(Z>1 .4 t Z<=7),RETURN(A39)) 

=IF(AND(AND(Z>7 t 2<=50)^AND(ROT<=3000 1 ROT>2000)),RETURN(lnterpolate!INTERPOLATE(3000,ROT l 2000,A34,A35 
=lF(AND(AND(Z>7 t 2<=50)^ND(ROT<=2000 ( ROT>1 000)), RETURN(lnterpolate!INTERPOLATE(2000,ROT l 1000, A35.A36 
=IF(AND<AND(Z>7Z<=50)AND(ROT<=1000,ROT>700)),RETURN{imerpolate!INTERPOLATE(1000,ROT,700,A36,A37))) 
=1F(AND(AND(Z>7 f Z<=5O) t AND(ROT<=7OO,ROT>5OO)).RETURN(lnterpolat0!INTERPOLATE(7OO 1 ROT,5OO,A37,A38))) 

=IF(AND(AND(Z>7 f Z<=50),AND(ROT<=500, ROT >0)),RETURN(" r/l Out of Range")) 

=IF(Z>50,RETURN("Z Out of Range")) 

=1F(Z<1 ,RETURN("Z Out of Range")) 

=IF(ROT>3000,RETURN("rA Out ol Range")) 

=IF(ROT<O.RETURN("rA Out of Range")) 

=RETURN("Macro ErO 


RgC92 

=ARGUMENT("Zb") 

= ARGUM ENT (" AOB92") 

=9.2456+0 32745*Zb-0.0013434*Zb A 2 +0.00000424 12*2^3 a/b = Infinity 

= 10.128+0.31493*Zb-0.00031213*Zb A 2-0.0000019266’Zb A 3 a/b = 3 

= 10.762+0.37273*Zb-0.00054319*Zb A 2-0.000002604*Zb A 3 &*> = 2 

= 11 .72+0.4 1376*Zb-0.0012646*Zb A 2+0.0000025354*Zb A 3 &*> = 1 5 

= 14, 84 1*0. 34501 *Zb+0. 000497 1 *Zb A 2 0 0000073339*Zb A 3 = 1 

=IF(Zfcx1.RETURN("Zb Out of Range. 1 ')) 


=IF(AND(AND(Zb>=1,Zb<=100).AOB92>3),RETURN(lnterpolate l iNTERPOLATE(1000000,AOB92,3.A59.A60))) 
=IF(AND(AND(Zb>=1 ,Zlx=100),AND(AOB92<=3,AOB92>2)).RETURN(!n1erpolale!INTERPOLATE(3 1 AOB92.2 1 A60 1 A61))) 
=IF(AND(AND(Zb>=1 ,Zlx=100),AND(AOB92<=2.AOB92>1 5)) RETURN(lnterpolale!INTERPOLATE(2.AOB92,1 .5.A61 ,A6 
=IF(AND(AND(Zb>=1.Zb<=100),AND(AOB92<=1.5.AOB92>=1i},RETURN(lnterpolate!INTERPOLATE(1.5,AOB92.1,A62,A 


3ruhnF:g 


=IF{Zb>100,RETUHNrZb Out of Range.")) 
=IF(AOB92<1 ,RETURN("a/b Out of Range.")) 
*RETURN(’Macro ErTI 


RgC9.4 

* ARGUM ENT (“Zb94") 

* ARGUM ENT ("AOB90 


*5.1455+0.171 18*Zb94 

a/b * Infinity i 

*2.5003*Zb94 A (0.4 1 62) 

a/b = Infinity i 

*5.7433+0.1 61 69*Zb©4 

a/b = 3 

*6.2273+0. 1 8037“ Zb94 

a/b = 2 

*6.9795+0. 19377 4 Zb94 

a/b = t .5 

*9.3546+0.25446“Zb94 

a/b = 1 


»IF(Zb94<1 .RETURNCZb Out of Range.*), IF(Zb94>1 00, RETURNfZb Out of Range."))) 

-IF(AND(AND(Zb94>=20^b94<=100),AOB94>3),RETURN(lnterpolateHNTERPOUTE(1000000 l AOB94 l 3,A78,A79))) 
*IF(AND<AND(Zb94>«lZb94<=20),AOB94>3),RETURN(lntefpolate!lNTERPOlATE(lOOOOOO f AOB94,3.A77,A79))) 
*IF(AND(AND(Zb94>«1.Zb94<*100) l AND{AOB94<=3,AOB94>2)) ) RETURN(tnterpolate!INTERPOLATE(3 I AOB94,2 t A79,A 
=IF(AND(AND(Zb94»*l Zb94<=100).AND(AOB94<=2.AOB94>1 .5)),RETURN(lnterpotaleHNTERPOLATE(2.AOB94.1 .5,Af 
*IF(AND(AND(Zb94>«l2b94<«lOO) 1 AND(AOB94<=1.5,AOB94>=1)).RETURN(lntefpolate!INTERPOUTE(1.5.AOB94 1 1,/ 

= RETURN (“Macro ErO 
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Interpolate 


INTERPOLATE 

=ARGUMENT("ROT1“) 

=ARGUMENT (“ROT2") 

=ARGUMENT('ROT3') 

=ARGUMENT('VAL1 ") 

=ARGUMENTfVAL3") 

=RETURN((VAL1-VAL3)*(ROT2-ROT3)/(ROT1-ROT3)+VAL3) 


Stress 


Stress Calculation Macros 

Equivalents -- Determines Pequivalent and Mequivalent. 
=ARGUMENT('Select') 

=ARGUMENT(’Radr) 

=ARGUMENT('Mlim") 

=ARGUMENT(’Mult’) 

= ARGU MENT (' Axlim") 

=ARGUMENT('Axult') 

Meqlim= Mlim+Axlim*Rad1/2 
MeqUlt= MUIt+AxUlt*Rad1/2 
Peqlim= Axlim+Mlim*2/Rad1 
Peqult= AxUlt+MUIt*2/Rad1 
=IF(Select1=1 ,R ETURN(Meqlim)) 
=IF(Select1=2,RETURN(MeqUlt)) 

=IF(Select1 =3,RETURN(Peqlim)) 

=IF(Select1 =4 ,R ETU RN(Peqult)) 

=RETURN(“Select not valid.") 


Stress_C -- Limit Stress. 

=ARGUMENT ("Mlim") 

=ARGUMENT ('Axlim') 

=ARGUMENT("Rad") 

=ARGUMENT ("Nost") 

=ARGUMENT ("Ast") 

=ARGUMENT('lst") 

= ARGUMENT (“bst") 

=ARGUMENT ("tsk") 

=ARGUMENT('tr) 

= ARGUM ENT ("bf ") 

^Mlim'RadyMOIlMoi^Rad.Nost.Ast.lst.tsk.bf.tl) 
=Axlim/(2*PI0*Rad*(bst*tsk+tl*bf)/bst+Ast*Nost) 
=RETURN(A31 +A32) 


ET_Stress -- Ultimate Stress. 
=ARGUMENT(*select") 
=ARGUMENT ("MUIt") 

= ARGUMENTS AxUlt") 

=ARGUMENT(“Rad") 

=ARGUMENT(’Ec") 

= ARGUMENT (“tsk") 
=ARGUMENT ("bst") 
=ARGUMENT ("Ast") 
=ARGUMENT("bt") 
=ARGUMENT("tl") 
=ARGUMENT(''d") 
=ARGUMENT (“nu") 
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Stress 


Nxult= -(MUtt/PI()/RacK'2+AxUlt/2/PI()/Rad) 


Nyuit=0 


kc= MAX(BruhnFig!FigC5.6(bst/tsk,1 ),BruhnFig!FigC9.1 (MIN(bst,d) A 2/Raci/tsk*( 

Fcr= -kc*PI0 A 2*Ec/(1 2‘(1 -nu*2))*(tsk/bst)*2 
AJd=bf*(tsk+tl) 

tsm=(Ast+Ald+(bst-bf)*tsk)/bst 

tb=(tsk+(tsk+tl))/2 

Fskp= 1Asm*(Nxult+nu*Nyult‘(tsm-tsk)/tsk) 

=SET.VALUE(B63,Fskp) 

Ald=bt*(tsk+tl) 

we= 1 ,7*tb*SQRT(ABS(Ec/Fskp)) 

=IF((Fskp)>=Fcr,RETURN(’No Skin Failure’)) 

Pb= 0.9*Fcr*tsk*(bst-we) 

Rot=Nxult*bst 
Pwe= Rot-Pb 


Nxp= Pwe/we 


tsm=(Ast+ Ald+(we-bf) *tsk)/we 

Fskp= 1 Asm*(Nxp+0.33*Nyult*(tsm-tsk)Ask) 

we=1 ,7*tb*SQRT(ABS(Ec/Fskp)) 

=IF(we>bst,SET.NAME('we",bst)) 

=IF(we<bf,SET.NAME(’Ald',we*(tsk+tl))) 

=IF(ABS(Fskp)-ABS(B63)>0.01,GOTO(A63)) 


Fst= 1Asm*(Nxp-0.33*Nyult) 

Fskin= 1Asm*(Nxp+0.33*Nyutt*(tsm-tsk)Ask) 
=IF(select=1,RETURN(Fskin),RETURN(we)) 
=RETURN(Fskp) 


Beam_Theory -- Ultimate Beam Theory. 
= ARGUM ENT ("select") 

=ARGUMENT ("Mult") 
=ARGUMENT("AxUlt") 
=ARGUMENT("Rad") 
=ARGUMENT("Ec") 

=ARGUMENT("tsk") 

= ARGU M E NT(” Nost “ ) 
=ARGUMENT("bst'') 


-31 503.5836277C 


Stress 


=ARGUMENT('Ast‘) 

=ARGUMENT('bf) 

=ARGUMENT(tT) 

=ARGUMENT(’d') 

=ARGUMENT ('nu') 

DA=2*PI()/Nost 

Ald= (tl+tsk)*bf 

tav= (tsk*(bst-bf)+Ald)/bst 

tb=(tsk+(tsk+tl))/2 

moist=0 

moisk= PIQ*RacK'3*tav 


=FORCCountM,Nost) 
d=Rad*SIN((Count-1 )*DA) 

=SET.NAME('moist*, moist + Ast*d A 2) 
=NEXTQ 


moi= moist+moisk 


Cent=0 

Area=2*PI0*Rad*tav+Nost*Ast 

Aeff=0 

Stress max= 0 
wemax=0 
Atrack= 0 
DC=0 


leff= 0 


StressA= AxU It/Area 

kc= MAX(BruhnFig!FigC5.6(bst/tsk,1),BruhnFig!FigC9.1(MIN(bst,d) A 2/Rad/tsk*( 
Fcr= -kc*PI() A 2‘Ec/(12*(1-nu A 2))*(tsk/bst) A 2 

=FOR(’Count’,1 ,Nost) 
d=Rad*S!N((Count-1 )*DA) 

Stress=- MUtt*(d-Cent)/moi-StressA 

we=1 ,7*tb*SQRT (ABS(Ec/Stress)) 

=IF(we>bst,SET.NAME("we“,bst)) 

=IF(Stress<Stressmax) 

=SET.NAME("Stressmax“, Stress) 

=SET.NAME(“wemax'',we) 

=END.IF() 

Ae=ABS((bst-we)*tsk*Fcr/Stress) 

Ald=bf*(tl+tsk) 
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Stress 


=IF(we<bf,SET.NAME("Ald",we*(tl+tsk))) 

tav=((we-bf)*tsk+Ald)/we 

=IF(we<bf,SET.NAME(*tav*,Ald/we)) 

=SET.NAME('Aeff’,Aeff+(we*tav+Ast+Ae)) 

=S ET. N AM E ('At rack* , Atrac k+(we*tav+ Ast + Ae) *d) 
=SET.NAME('leff\leff+(we*tav+Ast+Ae)*d A 2) 


=NEXTQ 


=SET.NAME(‘DC',Cent-Atrack/Aeff) 

=SET.NAMEfCent‘,Atrack/Aeff) 


=S ET. NAM E(* Area " ,Aeff) 

=SET.NAME('moi',leff+Aeff*DC A 2) 

=IF(ABS(DC)>0.1 ,GOTO(A122)) 

=1 F(select=1 , R ETU R N(Stressmax) , R ETU R N(wemax)) 
=RETURN(*Macro Error*) 
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APPENDIX D 
End-Fixity Coefficients 
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Fixity Coefficient, 
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Fixity Coefficient for a Column with Simply Supported 
Ends and an Intermediate Support of Spring Constant. - 







Fixity Coefficient, 




Where ( p) = Spring 
Constant Which Is Equal 
To The Number of Pounds 
Necessary To Deflect 
The Spring One Inch 
Extrapolated to Zero 


Deflection. 
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Fixity of a Column vitr. Two elastic, Symmetrically 
Placed Supports Having Spring Constants, p 
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APPENDIX E 


Bifurcation Buckling 


Source Code 


This program utilizes the rr.etkcds put forth in "Design 
Criteria for Axially Loaded Cylindrical Shells" and 
"Buckling of Bars, Plates, and Shells", by 
B . 0 . Almroth 


Real Ist,Ir # Jst, Jr,nu, L,C (6, , A<3, 3) ,AO(2,2) ,Al(3,3) ,Ny,p 

Real Lam(15),Nx(15,18) ,Nxp(15, 18) ,Ncl,Eta(15) ,Nwc,Ncr,mb 
Real Per, Pci (15, 15) , Pcla ( 15) 

Integer m,n 

CHARACTER TAB 
CHARACTER CR 
CHARACTER LF 
TAB = CHAR(9) 

CR = CHAR (13 ) 

LF = CHAR ( 14 ) 


pi = 3.14159 
p=17 , 3 

Stress = 19177,1 
E = 10.8E6 
nu = 0.33 
G = E/2/ (1+nu) 
bs = 10.832 
L = 40.08 

Open (unit =3, f ile= 1 Bifurcation. Dat 1 , s t at us= 'unknown* ) 
Write(3, '(/)*) 

Fpl = 5. 8*pi**2*E/12/ (l-nu**2) * (tsk/bs) **2 

tsk = 0.126 
Ast = 0.279 
1st = 0.04238 
Jst = 0.023682 
es = -0.77374 
R = 165 . 5-tsk/2 

twbar = tsk+Ast/bs 
tebar = (Ast+tsk*bsk) /bs 
Write(*,*)tebar 

Ar = 0.0 
Ir = 0.0 
dr = 40.08 
Jr = 0.0 
er = 0.0 


Cp = E*tsk/ (l-nu**2) 

D = E*tsk**3/12/(l-nu**2) 


*** Classical Bifuraction Buckling Analysis *** 

Do 10 1=1,6 
Do 11 J=1 , 6 
Cfi, j) = 0.000 
Continue 
Cont inue 


C(l, 1) = Cp+E*Ast/bs 



0 ( 1,2 

C(i,4; 


= nu*Cp 
= E*Ast ’es/bs 

0(2,1) = nu*Cp 
0(2,2) = Cp + E*Ar/dr 
0(2,5) = E*Ar*er/dr 

0(3,3) = G*tsk 

0(4,1) = E*Ast*es/bs 

0(4,4) = Cn-E/bs* (Ist+Ast*es**2) 

0(4,5) = nu*D 

0(5,2) = E*Ar*er/dr 
0(5,4) = nu*D 

0(5,5) = D + E/dr* { Ir+Ar*er* *2 ) 

C (6, 6) = 2* (1-nu) *D + G* (Jst/bs+Jr/dr) 

Write (3, 15) ((C(i,j), i=l, 6) , j=l, 6) 

15 Format (6f 12.1) 

Write(3,315) 

Do 20 i=l, 3 
Do 21 j=l, 3 
A( i, j ) = 0.0 
Al(i,j) = 0.0 
21 Continue 

20 Continue 

Do 30 i=l, 2 
Do 31 j=l,2 
Ao(i,j) = 0.0 
31 Continue 

30 Continue 


Do 101 m=l, 15 
Do 100 n=l, 15 


nn = n+3 

Lam(m) = (m*pi/L)**2 
Eta(n) = (n/R) **2 


Ad.l) = C(l, 1) *Lam(m) + C(3, 3) ’Eta(n) 
A{1,2) = (C(l,2) + C(3, 3) ) *m*pi*n/L/R 


A(l,3) = C(l,2)*m*pi/R/L+C(l,4)*Lar 
3 (C (1, 5) +2*C(3, 6) ) *m*pi*Eta (n) /L 


m) *m*pi,- C+ 


A ( 2 , 2 . = C(2,2: 

A(2,3 ; = (C(l, 5 
1C (2 , 5 ) ’Eta (n) *n. ?. 


Eta(n) + C(3, 3) m; 
+2*C(3, 6 ) ) ’Lamfiri' ’r. ?.+c (2, 2 


’n/R”2+ 


AO,;. = C (4 , 4 ) ’Lam (m) ”2+ <C 1 5 , 6 
1C (5, 5) ’Eta(n) **2+C(2, 2) /R**2+2*C!2 
32*C(1, 5) ’Lam(m) /R 


-2’C (4, 5) ’Eam(m) ’Eta (n)i- 
: P.’Eta in. - 


0RIGI*4AL PAG£ IS 
OF POOR QUALITY 
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c 

C211 

C 

C215 

C 

315 

C 

C415 


100 

101 


AO ( 1 » i ) 

= Add) 

Ao ( 1 , 2 ) 

= Add) 

Ao ( 2 , 1 ) 

= Add) 

Ao ( 2 , 2 ) 

= A(2, 2) 

Al (1,1) 

= Add) 

Al(l, 2} 

= Add) 

Aid, 3) 

= A(l, 3) 

Al (2, 1) 

= Add) 

Al(2,2) 

= A(2, 2) 

Al(2, 3} 

= A(2, 3) 

Al ( 3 , 1 } 

= A(l, 3) 

Al ( 3 , 2 ) 

= A(2, 3) 

Al ( 3 , 3 ) 

= A(3, 3) 

DETAo = 

Ao (1 ,1 ) * Ao ( 2 ,2 ) -Ao (2,1) *Ao (1,2) 


DETAl = Al(l, 1) MAI (2, 2) *Al (3, 3) -Al (3, 2) *A1 (2, 3) ) 
e-Al(l,2)*(Al(2, 1)*A1(3,3)-A1(3, 1)*A1(2,3) ) 

@+Al (1,3)* (Al (2,1) *Al (3,2) -Al (3,1) *A1 {2 , 2) } 

Ny = -p*R 

Nxp(m,n) = (l/Lam{m))*(DetAl/DetAo-Ny*Eta(n)> 

Nx (m, nn) =Nxp(m,nn-3) 


If (m .BQ. 1 .AND. n .EQ.l) Ncl=Nxp(m,n) 
If (Nxp (m,n) .LT. Ncl) Ncl = Nxp(m,n) 


Write {3, 211) m,n 
Format (2x, 12, 3x, 12, 3x f 12) 

Write(3,215) ( (Ao(i, j> , i=l,2) , j=l,2) 

Format (2fl2 .2) 

Write (3, 315) 

Format ( / ) 

Write (3, 415) ((AKi, j), i=l, 3) , j=l, 3) 

Format ( 3 f 2 0 .2) 

Write (3, • <2x, "NxpC '12, lh, , 12, ' •)= ",F20.4i'} m,n,Nxp(m,n) 

Continue 

Continue 


Write{*, ’ /* ‘Classical Bifurcation Buckling Allowable = 

(if 10. 2) ') Ncl 

Write (3, ‘Classical Bifurcation Buckling -llowarle = 

(if 10 . 2 ) ‘ ) Ncl 
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n C) 


**** wide Column Eucklir.r Analysis T ’’ r 


Cb44 = C(4,4)-C(l,4: "2/CU.l) 

Cb55 = C(5,5)-C(2,5 "2/C<2,2) 

Nwc = Cb44* (pi/L) **2 

Write{*, ' (' 'Wide Column Buckling Allowable 
@f 10.2) ’ ) Nwc 

Write(3, ' {' 'Wide Column Buckling Allowable 
@£10 .2, ) ' ) Nwc 


Det = C (1, 1) *C (2, 2) -C ( 1, 2) **2 

te = 1/sqrt (5 . 46* (Cb44+Cb55) *C (2 , 2) /Det ) 

ROte = R/te 

Phi = 6. 48/ROTe** (0.34371769) 

Write ( * , * ) ROTe , Phi 

Ncr = Nwc+phi* (Ncl-Nwc) 

Write ( * , ' ( ' 'Critical Buckling Allowable = 
@fl0.2, ) ■ ) Ncr 

Write{3, ' ( ' 'Critical Buckling Allowable = 
@f 10 . 2 , ) ' ) Ncr 


Close (unit=3) 

Stop 

End 


ORlGiWK PAOC K 
OF POOR QUALITY 


Output 


1805277.3 

503943.5 

.0 

-215235.7 

.0 

503943.5 

1527101.4 

.0 

.0 

.0 

.0 

.0 

511578.9 

.0 

.0 

-215235.7 

.0 

.0 

210811.6 

666.7 

.0 

.0 

.0 

666.7 

2020.4 

.0 

.0 

.0 

.0 

.0 


.0 

.0 

.0 

.0 

11584. CL- 



Nxp ( 

lr D = 

Nxp ( 

1, 2) = 

Nxp( 

1, 3 ) = 

Nxp{ 

1, 4) = 

Nxp( 

1, 5) = 

Nxp { 

1, 6) = 

Nxp( 

1, 7) = 

Nxp ( 

1, 8) = 

Nxp( 

1, 9) = 

Nxp ( 

1,10) = 

Nxp ( 

1,11) = 

Nxp ( 

1,12) = 

Nxp ( 

1,13) = 

Nxp ( 

1,14) = 

Nxp( 

1,15) = 

Nxp ( 

2 , 1) = 

Nxp{ 

2, 2) = 

Nxp( 

2, 3 ) = 

Nxp{ 

2, 4 ) = 

Nxp ( 

2, 5 ) = 

Nxp ( 

2, 6 ) = 

Nxp ( 

2 , 7 ) = 

Nxp( 

2, 8 ) = 

Nxp( 

2, 9 ) = 

Nxp{ 

2,10)= 

Nxp( 

2,11)= 

Nxp{ 

2,12)= 

Nxp ( 

2,13)= 

Nxp ( 

2,14)= 

Nxp ( 

2,15)= 

Nxp { 

3, 1) = 

Nxp ( 

3, 2) = 

Nxp ( 

3, 3 ) = 

Nxp ( 

3, 4 ) = 

Nxp { 

3, 5) = 

Nxp{ 

3, 6 ) = 

Nxp ( 

3, 7) = 

Nxp { 

3, 8 ) = 

Nxp ( 

3, 9) = 

Nxp ( 

3,10)= 

Nxp{ 

3,11)= 

Nxp ( 

3,12)= 

Nxp ( 

3,13)= 

Nxp ( 

3,14)= 

Nxp ( 

3,15)= 

Nxp ; 

4, 1) = 

Nxp ( 

4, 2 ) = 

Nxp [ 

4, 3 ) = 

Nxp t 

4, 4 ) = 

Nxp \ 

4, 5 ) = 

Nxp \ 

4 , 6) = 

Nxp [ 

4, 7 } = 

Nxp i 

4, 8 } = 

Nxp i 

4, 91 = 

Nxp { 

4,10}= 

Nxp { 

4,11)= 


9999.0313 

9685.4014 

9210.9961 

8635.0781 

8020.2241 

7421.9658 

6883.0952 

6432.3403 

6085.9292 

5850.4272 

5725.7290 

5707,5947 

5789.5723 

5964.3306 

6224.4766 

7329.5503 

7305.3750 
7266.5332 
7215.0884 
7153.7222 
7085.5566 
7013.9600 

6942.3750 
6874.1475 
6812.4028 
6759.9502 
6719.2285 
6692.2764 
6680.7324 
6685.8589 

11878.1318 
11871.6064 
11860.9795 
11846.6309 
11829.0645 
11808.8945 
11786.8398 
11763.6924 
11740.2871 
11717.5029 
11696.2227 
11677.3135 
11661. 624C 
11649.9580 
11643.0674 
19441.2734 
19433.136'? 
19432.9961 
19425.9902 
19417.2949 
19407.1348 
19395.79:2 
19383.5430 
19370.7320 
19357.7734 
19344.9833 


Qi mrmi pa&l ^ 

OF POOR QUALITY 
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Nxp{ 
Nxp ( 
Nxp ( 
Nxp { 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp { 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp { 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp ( 
Nxp { 
Nxp { 
Nxp ( 
Nxp ( 
Nxp { 
Nxp ( 
Nxp ( 
Nxp { 


' - n ' _ 

19332. 7S52 

4,13) = 

19321. 5S01 

4, 14) = 

19311.7734 

4,15) = 

19303.7695 

5, 1) = 

29494.0020 

5, 2 ) = 

29491.9316 

5, 3 ) = 

29488.5430 

5, 4 ) = 

29483.8848 

5, 5) = 

29478.0527 

5 , 6) = 

29471.1543 

5, 7 ) = 

29463.3281 

5, 8) = 

29454.7402 

5, 9) = 

29445.5605 

5,10)= 

29435.9785 

5,11)= 

29426.1934 

5,12)= 

29416.4414 

5,13)= 

29406.9414 

5,14)= 

29397.9355 

5,15)= 

29389.6680 

6, 1) = 

41906.3320 

6, 2 ) = 

41904.7266 

6, 3) = 

41902.0664 

6, 4) = 

41898.3984 

6, 5) = 

41893.7773 

6, 6) = 

41888.2695 

6, 7 ) = 

41881.9609 

6, 8) = 

41874.9375 

6, 9) = 

41867.3008 

6,10)= 

41859.1953 

6,11)= 

41850.7305 

6,12)= 

41842.0508 

6,13)= 

41833.3008 

6,14)= 

41824.6367 

6,15)= 

41816.2188 

7, 1) = 

56633.6836 

7, 2 ) = 

56632.3086 

7, 3) = 

56630.0234 

7, 4 ) = 

56626.8633 

7, 5) = 

56622.8477 

7, 6) = 

56618.0430 

7, 7) = 

56612.4961 

7-, 8 ) = 

56606.2617 

7, 9) = 

56599.4336 

7,10)= 

56592.0703 

7,11)= 

56584.2813 

7,12)= 

56576.1328 

7,13)= 

56567.7500 

7,14)= 

56559.2500 

7,15)= 

56550.7305 

8, 1> = 

73657.4219 

8, 2 ) = 

73656.1563 

8, 3 ) = 

73654.0859 

8, 4 ) = 

73651.2109 

3, 5) = 

73647.5469 

3, 6 ) = 

73643.1484 

3, 7) = 

73638.0313 

3, 8) = 

73632.2578 

3, 9) = 

73625.8828 

8,10)= 

73618.9531 

3,11)= 

73611.5391 

3,12)= 

73603.7109 

5,13)= 

73595.5234 

3,14)= 

73587.0938 

3,15)= 

73578.4844 

?, 1) = 

92968.5938 

9, 2 ) = 

92967.4297 


QrtlSmi FAOE !5 
OF POOR QUALITY 


Nxp( 9, 

3 ) = 

Nxp( 9, 

4) = 

Nxp( 9, 

5) = 

Nxp( 9, 

6) = 

Nxp( 9, 

7) = 

Nxp( 9, 

8) = 

Nxp ( 9, 

9) = 

Nxp ( 9, 

10) = 

Nxp( 9, 

11) = 

Nxp( 9, 

12) = 

Nxp( 9, 

13) = 

Nxp( 9, 

14) = 

Nxp ( 9, 

15) = 

Nxp(10, 

1> = 

Nxp{10, 

2) = 

Nxp{10, 

3) = 

Nxp{10, 

4) = 

Nxp(10, 

5) = 

NxpdO, 

6) = 

NxpdO, 

7) = 

Nxp(10, 

8) = 

Nxp(10, 

9) = 

Nxp(10, 

10) = 

NxpdO, 

11) = 

NxpdO, 

12) = 

NxpdO, 

13) = 

NxpdO, 

14) = 

NxpdO, 

15) = 

Nxp (11, 

1) = 

Nxpdl, 

2) = 

Nxp (11, 

3) = 

Nxpdl, 

4) = 

Nxpdl, 

5) = 

Nxpdl, 

6) = 

Nxp(ll, 

7) = 

Nxp(ll, 

8) = 

Nxp(ll, 

9) = 

Nxp(ll, 

10) = 

Nxpdl, 

1D = 

Nxpdl, 

12) = 

Nxpdl, 

13) = 

Nxpdl, 

14) = 

Nxpdl, 

15) = 

Nxp (12, 

1> = 

Nxp ( 12 , 

2) = 

Nxp (12, 

3) = 

Nxp ( 12 , 

4 ) = 

Nxp (12, 

5 ) = 

Nxp (12, 

6) = 

Nxp (12, 

7 ) = 

Nxp (12, 

8 ) = 

Nxp ( 12 , 

9 ) = 

Nxp ( 12 , 

10) = 

Nxp (12, 

11) = 

Nxp ( 12 , 

12) = 

Nxp ( 12 , 

13 ) = 

Nxp (12 , 

14) = 

Nxp (12, 

15 ) = 

Nxp (13, 

1 ) = 

Nxp (13, 

2 ! = 

Nxp (13, 

3 ) = 

Nxp (13, 

4 ) = 

Nxp (13, 

5 ) = 

Nxp (13, 

6) = 

Nxp (13, 

7 ) = 

Nxp (13, 

8 ) = 


92965 . 4844 
92962.7813 
92959.3516 
92955.1797 
92950.3438 
92944.8672 

92938.7578 
92932.0859 
92924.9141 
92917.2734 
92909.2109 
92900.8125 
92892.1484 

114562.5625 
114561.4453 
114559.5938 
114557.0234 
114553.7031 

114549.7344 
114545.0547 

114539.7578 

114533.8516 

114527.3594 
114520.3359 
114512.8203 

114504.8516 
114496.4531 

114487.7578 

138436.6719 

138435.5781 
138433.7969 

138431.2969 
138428.0781 
138424,2031 

138419.6719 
138414.5313 

138408.7344 
138402.3438 
138395.4688 

138388.0000 

138380.0938 
138371.7500 
138362.9688 
164589.2188 
164588.1563 

164586.3906 
164583.9844 
164580.8281 
164577.0469 
164572.6250 
164567.5469 
164561.8906 

164555.5625 
164548.7188 

164541.3594 
164533.4844 

164525.0938 

164516.3594 
193019.3281 

193018.2969 

193016.5781 
193014.1875 
193011.1406 

193007.3906 

193003.0000 
192998.0313 



:ixp{13, 9) = 
Nxp (13,10) = 
Nxp(i3, ii)= 
Nxp (13,12) = 
Nxp( 13 , 13) = 
Nxp (13,14) = 
Nxp (13, 15) = 
Nxp (14, 1) = 
Nxp (14, 2) = 
Nxp (14, 3 ) = 
Nxp (14, 4 ) = 
Nxp (14, 5)= 
Nxp { 14 , 6)= 
Nxp { 14 , 7 ) = 
Nxp (14, 8 ) = 
Nxp (14, 9 ) = 
Nxp (14, 10) = 
Nxp (14,11)= 
Nxp(14, 12)= 
Nxp (14, 13)= 
Nxp ( 14 , 14 ) = 
Nxp (14 , 15) = 
Nxp { 15 , 1)= 
Nxp (15, 2 ) = 
Nxp (15, 3 ) = 
Nxp (15, 4 ) = 
Nxp (15, 5)= 
Nxp (15, 6)= 
Nxp(15, 7 ) = 
Nxp ( 15 , 8 ) = 
Nxp (15, 9) = 
Nxp(l5,10)= 
Nxp( 15, 11 ) = 
Nxp (15, 12) = 
Nxp (15,13)= 
Nxp (15,14) = 
Nxp (15, 15) = 


192992.3906 

192986.2031 

192979.3750 

192972.0625 
192964.2188 
192955.8594 

192947.0625 
223726.2344 
223725.1719 
223723.5156 
223721.0938 

223718.1406 

223714.4531 
223710.1250 
223705.1875 
223699.6094 

223693.4531 
223686.6875 

223679.3750 
223671.5781 

223663.2031 

223654.3750 
256709.4375 
256708.4219 
256706.7656 
256704.4688 
256701.4063 

256697.8125 
256693.5000 
256688.5938 

256683.1406 
256676.9688 
256670.2656 
256662.9844 
256655.1563 

256646.8125 
256637.9375 


Classical Bifurcation Buckling Allowable = 5707.59 

Wide Column Buckling Allowable = 1137.54 

Critical Buckling Allowable = 3231.44 


CL - 2 _ 
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APPENDIX F 


Margin of Safety Calculations 


93 


IBRffflONftLV HU* 


PRECEDING PAGE BLANK NOT FILMED 



General Instability 
Shanley 


Margins of Safety 


M C r = 272,794,759 in-lb 

P eq = [Bruhn, “Missiles,” page E 1 .98] 

P eq = 3,296,613 lb 

P. a 

Load per stringer column = ^~ 34,339.72 lb 
Area per column = 0.8183 in 2 

Stress (allowable) = -^ = ^0 8^83^ = 41,964.7 Ib/in 2 

& ALLOWABLE'*' ® YIELD PcR = Pcy~ 38,000 lb/in 2 

xjo 38,000 1 non/:'* 

MS = gA/t - 1 = 0.2062 . 

31,504 
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Margins of Safety 


Local Stringer Buckling 


F C r = 51,000 lb/in 2 


Limit Stress = 13,730 lb/in 2 


MS = 


51,000 

13,730 


-1=2.715. 


Column Buckling and Stringer Crippling 

Lowest critical stress in the column results from Euler buckling. 

F CR =15,917 lb/in 2 

c *Eiiler 

MS = _ 1 = _0 4947 

31,503 



Margin of Safety 


Skin Buckling: Combined Compression and Shear 


MS = — 77777 [Bruhn, “Aircraft,” page C5. 1 1 ] 

Rc^RlURf\ m 

n fc p Js 

Kc ~fT ’ Ks ~fT 


_ 13,729.9 
c " 7,823 


1.7551 


*s = 


713.1 

8,160 


0.0874 


MS = 


1.7551+[(1.7551) 2 +4(0.0874) 2 ] 1/2 


MS = -0.4316 
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